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I.  INTRODUCTION 

A  jet  sheet  directed  downwards  from  the  trailing  edge  of  a  wing,  as 
shown  in  Figure  1,  significantly  increases  the  wing  lift,  especially  at 
low  speeds.  Lift  is  generated  both  as  a  direct  reaction  to  the  jet 
impulse  and  as  a  result  of  changes  in  the  wing  surface  pressure  distribu¬ 
tion  induced  by  deflection  of  the  main  stream.  This  phenomenon  is  called 
the  jet  flap  effect  by  analogy  to  the  action  of  a  mechanical  flap.  It 
makes  some  contribution  to  the  forces  induced  by  all  powered  lift  systems 

In  addition  to  generating  a  total  lift  on  the  wing  greater  than  the 
vertical  component  of  the  jet  impulse,  the  jet  flap  also  induces  a 
thrust  greater  than  the  horizontal  component  of  the  jet  impulse.  This 
was  discovered  during  testing  of  jet  flapped  airfoils  at  the  British 


Figure  1.  Jet  Flap  Thrust  Recovery  Hypothesis 
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National  Gas  Turbine  Establishment.  As  a  result  of  this  discovery, 
Stratford  (ref  1)  suggested  the  hypothesis  that,  in  ideal  flow,  the 
total  jet  thrust  is  always  recovered  as  a  horizontal  force,  regardless 
of  the  initial  deflection  angle,  because  the  jet  is  ultimately  turned 
in  the  direction  of  the  free  stream.  He  did  not  attempt  a  formal  proof 
of  this  hypothesis,  but  reasoned  from  first  principles  that  the  force 
on  the  airfoil  must  be  equal  but  opposite  to  the  momentum  of  the  deflected 
jet. 

This  thrust  recovery  hypothesis  was  received  with  considerable 
skepticism,  but  the  existence  of  the  phenomenon  was  established  by  tests 
performed  at  Stanford  University,  when  Foley  (ref  2)  measured  almost 
complete  thrust  recovery  for  jet  deflection  angles,  6,  up  to  60°.  How¬ 
ever,  Stratford  (ref  3)  and  Williams,  et  al  (ref  4)  measured  only  partial 
recovery  at  all  deflection  angles,  with  large  decreases  in  recovery  as 
the  jet  deflection  angle  increased  to  90°.  In  subsequent  tests  with 
Foley's  model,  Tsongas  (ref  5)  also  found  a  decrease  in  recovery  when 
the  deflection  angle  was  increased  beyond  60°,  although  Hynes  (ref  6) 
found  that  the  flow  separated  from  the  tunnel  wall  in  this  case. 

The  loss  in  thrust  recovery  has  been  attributed  to  various  three- 
dimensional  and  real  fluid  effects:  Maskell  and  Spence  (ref  7)  pointed 
out  that  deflection  of  the  jet  wake  by  the  downwash  behind  a  finite  wing 
would  reduce  the  recovery,  even  in  ideal  flow.  Stratford  (ref  3)  and 
Wygnanski  (ref  8)  reasoned  that  the  thrust  recovery  is  decreased  in  real 
fluids  by  the  jet  entrainment,  which  induces  a  form  drag  on  the  airfoil; 
they  suggested  that  the  entrainment  and,  therefore,  the  jet  drag  increases 
with  the  jet  deflection  angle.  More  recently,  Wilson  (ref  9)  used  an 
empirical  entrainment  function  to  compute  the  jet  drag  and  show  that  it 
is  on  the  order  of  the  observed  loss  in  recovery.  On  the  other  hand, 
Williams,  et  al  (ref  4)  argued  that  flow  separation  on  the  airfoil  would 
also  cause  a  loss  of  thrust,  while  Tsongas  (ref  5)  suggested  that  separa¬ 
tion  of  the  main  stream  from  the  jet  flap  was  the  loss  mechanism.  Thus, 
while  there  is  clear  evidence  of  a  thrust  recovery  phenomenon,  the  real 
fluid  effects  which  reduce  the  recovery  are  incompletely  understood. 


Various  analytic  methods  have  been  developed  to  predict  the  increments 
of  aerodynamic  lift  and  pitching  moment  induced  by  the  jet  flap.  These 
methods  are  based  on  Spence's  (ref  10)  now  classical  jet  flap  theory,  in 
which  the  inertia  of  the  jet  sheet  is  related  to  the  strength  of  an 
equivalent  vortex  sheet.  The  original  theory  has  been  extended  to 
include,  for  example,  three  dimensional  effects  (Maskell  and  Spence, 
ref  7),  ground  effects,  (Lissaman,  ref  11),  and  jet  oscillations  (Simmons 
and  Platzer,  ref  12).  In  addition,  the  theory  has  been  applied  to  over¬ 
wing-blowing  (Lan,  ref  13)  and  ejector  wing  (Wilson,  ref  9)  configurations. 

Although  good  predictions  of  the  lift  and  moment  have  been  obtained, 
none  of  these  methods  predict  the  loss  of  thrust  recovery.  In  current 
practice,  complete  recovery  is  assumed  for  calculating  lift,  while  an 
empirical  function  derived  from  the  Stanford  data  is  used  to  specify  a 
reduction  in  the  thrust  recovery.  Entrainment  functions  have  been  use¬ 
ful  for  estimating  the  magnitude  of  the  thrust  loss,  but  such  functions 
treat  just  half  of  the  jet-/airfoil  interaction,  the  jet  drag  induced  on 
the  airfoil;  the  other  half  of  the  interaction,  the  equal  but  opposite 
reduction  of  the  jet  thrust,  is  neglected.  Since  the  reduction  in  thrust 
affects  the  computed  position  of  the  jet,  which  in  turn  affects  the  com¬ 
puted  force  on  the  airfoil,  a  method  for  predicting  the  effect  of  the 
jet  flap  should  include  this  interaction. 

The  jet  flap  methods  are  further  limited  by  linearizing  assumptions 
to  low  values  of  the  jet  thrust  coefficient  and  small  airfoil  and  jet 
deflection  angles.  Thus,  they  cannot  be  used  to  predict  the  takeoff  or 
landing  performance  of  V/STOL  aircraft  for  which  current  interest  in  the 
jet  flap  effect  is  greatest.  Addessio  and  Skiftad  (ref  14)  recently 
addressed  some  of  these  low  speed  effects,  but  made  the  usual  assumption 
of  complete  thrust  recovery. 

Such  semi-empirical  methods  are  only  useful  for  small  variations  from 
the  original  data  base.  An  essentially  complete  range  of  jet  deflection 
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angles  has  been  tested,  but  all  of  the  data  are  for  jet  thrust  co- 
effocients  less  than  unity.  Since  the  thrust  coefficient  is  much 
greater  than  unity  for  most  of  the  V/STOL  flight  regime,  there  are  no 
data  for  either  the  force  on  the  airfoil  or  the  trajectory  of  the  jet 
sheet  in  this  case. 

The  purpose  of  this  report  is  to  describe  progress  in  developing  a 
theory  of  jet  flap  thrust  recovery  at  the  low  speeds  and  large  deflec¬ 
tion  angles  characteristic  of  V/STOL  lift  systems.  A  combination  of 
analysis  and  testing  has  been  utilized  to  develop  this  theory.  A  viscid/ 
inviscid  interaction  analysis  of  the  jet  flap  has  been  devised  to 
compute  the  contribution  of  the  jet  drag  to  the  reduction  in  thrust 
recovery.  In  addition,  wake  survey  and  airfoil  surface  pressure  measure¬ 
ments  have  been  made  for  a  two-dimensional,  jet  flapped  airfoil  over  a 
range  of  jet  thrust  coefficients.  In  the  next  section  of  this  report, 
the  mechanism  of  thrust  recovery  and  the  way  in  which  entrainment  and 
separation  reduce  the  recovery  are  discussed.  The  jet  drag  analysis  is 
then  presented  in  the  following  section.  The  experimental  apparatus  and 
test  procedure  are  described  in  the  section  after  that.  In  the  last 
section,  the  analytical  results  are  compared  to  the  measured  thrust 
recovery.  It  is  concluded  that  the  thrust  recovery  is  nearly  complete 
for  small  values  of  the  jet  deflection  angle;  but  for  larger  angles, 
the  recovery  decreases  as  the  thrust  coefficient  increases.  For  values 
of  the  thrust  coefficient  on  the  order  of  unity,  the  thrust  loss  is  due 
to  jet  drag;  for  larger  values,  the  loss  is  increased  by  flow  separation 
from  the  airfoil. 


II.  JET  FLAP  THEORY 


MECHANISM  OF  THRUST  RECOVERY 

The  thrust  recovery  is  generated  by  equal  but  opposite  pressure 
forces  induced  on  the  airfoil  and  jet  sheet.  According  to  Spence’s 
(ref  10)  theory,  the  inertial  forces  which  resist  deflection  of  the  jet 
sheet  are  balanced  by  a  pressure  difference  in  the  external  flow  along 
the  upper  and  lower  boundaries  of  the  jet.  This  pressure  difference 
can  be  related  to  the  strength  of  an  equivalent  vortex  sheet  located 
along  the  jet  axis.  Then,  in  mathematical  terms,  the  vorticity  in  a 
section  of  the  jet  sheet  induces  an  upwash  velocity  of  magnitude, 
rj/2»T,  on  a  vortex  segment  bound  in  the  wing,  as  shown  in  Figure  2. 

An  upstream  force,  p rw(rj/2irr ) ,  is  therefore  induced  on  the  wing. 
Similarly,  the  wing  vortex  inducts  a  downwash,  rw/2irr  on  the  jet  vortex 
which  then  experiences  an  equal  force  prj(rw/2irr)  directed  downstream. 


Fieure  2.  Mutually  Induced  Forces  on  a  Vortex  Pair 


It  is  the  integrated  effect  of  all  such  vortex  pairs  which  results 
in  the  apparent  "thrust  recovery."  Figure  3  is  a  sketch  of  the  forces 
on  the  airfoil,  jet,  and  a  surrounding  control  volume.  As  seen  in  the 
figure,  the  jet  is  not  actually  deflected,  as  by  a  vane.  Rather,  the 
interaction  between  the  airfoil  and  jet  induces  a  horizontal  force,  pwy 
on  the  jet  and  a  corresponding  reaction  force  on  the  airfoil.  Simultan¬ 
eously,  the  interaction  of  the  main  stream  with  the  jet  cancels  the  verti¬ 
cal  component  of  the  jet  thrust  and  produces  an  equal  reaction  on  the  main 
stream.  Thus,  the  horizontal  thrust  of  the  jet  is  increased,  without  any 
reduction  in  the  vertical  reaction  force  on  the  airfoil. 

If  the  entrainment  of  the  jet  is  neglected,  Spence's  theory  predicts 
recovery  of  the  total  jet  thrust  as  a  horizontal  force.  This  can  be  shown 
from  the  equilibrium  of  the  forces  which  then  act  on  the  jet.  Consider 
the  segment  of  the  jet  shown  in  Figure  4.  The  forces  due  to  the  main 
stream  and  the  wing  vortex  system  deflect  the  jet  through  an  angle,  9, 


T  sin  9 


T  cos  6 


Spw7  +  T  cos  6 


L  +  T  sin  9 


Figure  3.  Forces  Induced  on  an  Airfoil  and  Jet  Flap 


Figure  4.  Forces  on  a  Section  of  the  Jet  Sheet 


as  shown.  The  initial  thrust  of  the  jet  is  T^  and  the  final  thrust  is 
T2*  Since  entrainment  has  been  neglected,  the  jet  will  be  a  streamline 
of  the  flow  and  the  pressure  force,  P,  must  therefore  act  normal  to  the 
segment . 

The  summation  of  forces  in  the  direction  of  T2  yields 

t2  =  T1  cos  9  +  P  sin  2-1 

The  summation  of  forces  normal  to  the  direction  of  T2  gives 

P  cos  j  «  Tj  sin  e  2-2 
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Solving  this  equation  for  P  and  substituting  in  Equation  2-1  gives 
T2  ■  Ti*  Therefore,  the  thrust  of  the  jet  is  conserved.  Because  the 
jet  consists  of  a  series  of  these  segments,  its  total  thrust  is 
ultimately  recovered  in  the  far  wake  as  a  horizontal  force. 

The  horizontal  force  induced  on  the  jet  is  the  integral  of  the 
corresponding  pressure  difference  across  the  jet;  that  is, 

oo 

Ipwy  -  f  dP  sin  6  d0  2-3 

Jo 

According  to  Spence  (ref  10) ,  dP  =  T/R.  Since  T  is  constant  and 
ds  =  Rd9,  integration  yields 

IpwY  =  T(1  -  cos  9)  2-4 

The  sum  of  this  force  and  the  horizontal  component  of  the  initial  jet 
momentum  gives  the  total  jet  thrust.  The  vertical  force  on  the  jet 
is  similarly  found  to  be  T  sin  8. 

This  thrust  recovery  mechanism  was  not  verified  by  the  Stanford 
experiments:  Quambeck  (ref  15)  integrated  the  surface  pressure  distribu 
tions  as  the  jet  deflection  angle  was  varied  and,  finding  that  the 
horizontal  force  did  not  change,  he  concluded  that  the  entire  jet  thrust 
was  felt  as  an  internal  pressure  on  the  nozzle.  However,  in  these  tests 
the  jet  was  deflected  by  utilizing  the  Coanda  effect  (Metral  and  Zerner, 
ref  16)  to  turn  it  over  a  short  flap.  The  reaction  to  the  force  which 
deflects  the  jet  over  the  flap  appears  as  a  reduction  in  the  static 
pressures  on  the  Coanda  surface.  The  resulting  pressure  drag  is  equal 
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but  opposite  to  the  thrust  induced  by  the  jet  flap  because  the  force 
initially  required  to  deflect  the  jet  over  the  flap  must  be  the  same 
as  the  force  necessary  to  bend  it  back  in  the  direction  of  the  free 
stream.  As  a  result,  the  net  horizontal  force  is  independent  of  the 
jet  deflection  angle.  Today,  it  is  no  longer  doubted  that  the  thrust 
recovery  is  induced  on  the  airfoil  surface,  primarily  as  a  leading  edge 
suction. 

thrust  loss  hypotheses 


Real  fluid  effects  reduce  the  thrust  recovery  through  two  separate 
mechanisms.  The  very  large  suction  peak  generated  near  the  leading  edge 
of  the  airfoil  increases  the  possibility  of  flow  separation  in  this 
region.  If  the  flow  does  separate,  the  resulting  loss  of  suction 
reduces  the  thrust  on  the  airfoil,  while  entrainment  of  the  low  momentum 
wake  fluid  by  the  jet  causes  a  corresponding  loss  of  thrust.  Separation 
becomes  more  likely  as  the  jet  thrust  coefficient  and  deflection  angle 
are  increased. 

However,  even  if  the  flow  does  not  separate  from  the  airfoil, 
entrainment  reduces  the  thrust  recovery.  A  jet  drag  is  induced  by  the 
mutuaLLy  induced  forces  on  the  vortex  sheet  bound  in  the  wing  and  a  sink 
distribution  which  represents  the  entrainment  of  the  jet.  The  origin  of 
the  force  on  the  wing  is  relatively  easy  to  understand;  the  jet  sinks 
change  the  surface  pressure  distribution  so  as  to  increase  the  wing 
drag.  There  are  two  contributions  to  this  force,  as  shown  in  Figure  5. 
The  jet  drag  due  to  thickness  originates  in  the  interaction  between  the 
jet  sinks  and  the  source  distribution  which  defines  the  airfoil  thick¬ 
ness.  The  jet  drag  due  to  lift  originates  in  the  interaction  between 
the  jet  sinks  and  vortex  distribution  which  determines  the  wing  lift. 

The  origin  of  the  reaction  force  on  the  jet  may  not  be  as  apparent. 
This  force  is  conceptually  similar  to  the  drag  experienced  by  a  sink 
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SINK-SINK  DRAG  DUE  TO  MIXING 


VORTEX-SINK  DRAG  DUE  TO  MIXING 


Figure  5.  Origin  of  the  Mixing  Loss 

in  an  external  stream.  However,  the  sink/vortex  interaction  is  only 
an  irrotational  simulation  of  the  jet  mechan'sm.  The  fluid  entrained 
into  the  jet  becomes  rotational,  so  that  au  understanding  of  how  the 
reaction  force  actually  develops  requires  consideration  of  the  jet 
mixing  process. 

This  process  is  basically  an  inelastic  collision  between  the  jet  and 
surrounding  fluid.  As  such,  jet  mixing  is  governed  by  the  same  laws  of 
momentum  and  energy  conservation  as  simple  collisions  between  discrete 
particles.  If  the  mixing  occurs  in  a  region  of  constant  static  pres¬ 
sure,  the  thrust  of  the  jet  is  conserved.  But  if  the  jet  passes  through 
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a  region  where  the  static  pressure,  P*,  is  different  from  the  undis¬ 
tributed  pressure  at  infinity,  the  mixing  changes  the  momentum  flux. 
The  following  simple  analysis  illustrates  this  phenomenon. 


Inside  the  region  where  the  pressure  is  P*,  the  velocities  in  the 
jet  and  coflowing  stream,  before  mixing,  are  assumed  to  be 


U»*  = 


2-5 


and 


Us*  = 


2-6 


in  which  PQ  is  the  stagnation  pressure  of  the  jet  and  U,*,  is  the  velocity 
of  the  coflowing  stream  at  infinity.  Momentum  is  conserved  during  the 
mixing  process  itself,  so  that 


thjUj*  +  msUs*  =  (mj  +  ms)  Um*  2-7 

in  which  mj  is  the  initial  mass  flux  of  the  jet,  ms  is  the  quantity  of 
mass  entrained  by  the  jet,  and  Um*  is  the  average  velocity  of  the  mixed 
flow. 


The  velocity  of  the  mixed  stream  changes  as  it  passes  into  the  region 
of  undisturbed  pressure;  the  velocity  becomes 


um 


2-8 


The  ratio  of  the  final  thrust  of  the  mixed  stream,  (iftp  +  &s  )  Urai  to  the 
thrust,  ifipUj ,  obtained  by  an  isentropic  expansion  of  the  jet,  without 
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mixing,  to  the  final  pressure  may  be  evaluated  by  substituting  in  turn 
for  Um,  then  U{£,  and  finally  for  Uj*  and  Us*.  Performing  these  sub¬ 
stitutions  yields  for  this  ratio 

4>  =  £l+2M  <  (  1+h)1/2  (U2+H)V2  -  H>  +M2U2J  1/2  -  MU  2-9 

in  which  M  =  ms/mp,  is  the  entrainment  ratio,  U  =  U»/Uj  is  the  velocity 
ratio,  and  H  =  2(PaB  -  P*)/  Uj2  is  the  normalized  pressure  change. 

If  there  is  no  mixing  between  the  jet  and  the  coflowing  stream,  then 
M  =  0  and  the  thrust  of  the  jet  is  conserved;  that  is,  <f>  =  1.  Similarly, 
if  there  is  no  pressure  change,  then  H  =  0  and  the  solution  reduces  to 
that  for  free  jet  mixing;  the  thrust  is  conserved  in  this  case,  also. 
However,  if  the  pressure  in  the  mixing  region  is  less  than  the  pressure 
at  infinity,  then  H  <  0  and  the  thrust  of  the  jet  increases.  For  the 
case  in  which  =  0,  the  thrust  ratio  has  the  simple  form, 

<t>  -  (1  +  M)V2  2-10 

and  the  jet  thrust  increases  with  increasing  entrainment.  This  is  the 
basic  mechanism  of  ejector  thrust  augmentation,  which  has  been  discussed 
in  greater  detail  by  Bevilaqua  (ref  17). 

If  the  pressure  in  the  mixing  region  is  above  ambient,  the  thrust 
of  the  jet  is  reduced.  For  the  case  in  which  Us*  =  0,  the  thrust  ratio 
also  has  a  simple  form 

^  a  1  -  MU  2-11 

and  in  this  case  the  jet  thrust  decreases  with  increasing  entrainment. 
Because  there  is  a  region  of  increased  static  pressure  behind  the 
trailing  edge  of  the  wing,  the  mixing  loss  mechanism  will  reduce  the 


jet  thrust.  According  to  Newton's  law  of  action  and  reaction,  the  net 
loss  of  thrust  must  be  equal  to  the  jet  drag  on  the  wing. 

The  thrust  loss  is  specified  by  a  thrust  recovery  coefficient,  r, 
which  is  defined  to  be  a  fraction  of  the  total  thrust;  that  is, 

rT  =  S  +  T  cos  60  2-12 

in  which  T  cos  0o  is  the  horizontal  cnrsp<*n *  of  the  jet  thrust  within 

the  nozzle,  and  S  is  the  horizontal  i.»t"  the  airfoil  surface 

pressure  distribution.  In  ideal,  t-wt-  'tfitsa -• -ooal  flow  S  =  T(1  -  cos  dQ ) 
and  so  r  =  1.  If  a  Coanda  flap  is  u.?e4  f^Jlect  the  jet,  the  pressure 
drag  due  to  the  Coanda  effect,  T( )  -  S*Q/,  is  included  in  S. 

Therefore,  the  recovery  is  defined 

rT  =  S  +  T  2-13 

In  this  case,  S  =  0  for  ideal,  two-dimensional  flow  so  that  r  =  1, 
again.  The  evaluation  of  the  recovery  coefficient  will  be  discussed  in 
greater  detail  in  Section  4.1,  on  data  reduction. 

In  order  to  determine  the  relative  importance  of  separation  and  jet 
drag  in  reducing  the  thrust  recovery,  the  jet  drag  will  be  computed  for 
unseparated  flow  and  compared  to  the  drag  actually  obtained  on  a  wind 
tunnel  model  under  the  same  conditions.  Differences  between  these  cases 
will  be  studied  to  quantify  the  loss  mechanisms.  In  the  next  section, 
the  method  of  calculating  the  jet  drag  will  be  developed  from  the 
principles  outlined  in  this  section. 


III.  COMPUTATIONAL  METHODS 


The  Navier-Stokes  equations  for  steady  fluid  motion  are  elliptic, 
which  means  that  the  domain  of  influence  of  a  point  disturbance  is  the 
entire  flow  volume.  That  is,  pressure  and  stress  gradients  transmit 
the  effect  of  local  disturbances  to  every  other  point  in  the  flow. 

Thus,  the  forces  induced  on  the  airfoil  depend  on  the  jet  position, 
while  the  jet  position,  in  turn,  depends  upon  the  forces  on  the  airfoil. 
The  basis  of  our  analysis  is  the  division  of  the  flow  field  into  two 
regions:  the  region  of  turbulent  mixing  within  the  jet,  and  the  inviscid 
region  surrounding  it.  The  equal  but  opposite  forces  on  the  airfoil 
and  jet  will  be  computed  by  matching  a  viscous  solution  for  the  jet  to 
a  potential  flow  solution  for  the  airfoil  and  jet  flap.  The  significant 
elliptic  effects  are  then  transmitted  through  the  pressure  field  of  the 
external  flow. 

For  the  viscous  jet  solution,  a  finite  difference  analysis  utilizing 
a  two-equation  turbulence  model  was  developed.  The  basic  scheme 
was  that  devised  by  Patankar  and  Spalding  (ref  18)  for  thin  shear 
layers:  the  full  Navier-Stokes  equations  are  reduced  to  a  simpler  set 
by  the  assumption  that  there  is  a  primary  direction  of  flow  (along  the 
jet),  and  that  the  diffusion  of  fluid  properties  is  negligible  in  that 
direction.  The  effects  of  curvature  on  the  development  of  the  jet  will 
be  treated  as  a  perturbation  of  this  basic  scheme.  A  jet  flap  panel 
method,  in  which  the  thrust  of  the  jet  is  not  assumed  to  be  constant, 
was  used  for  the  inviscid  solution. 

INVISCID  JET  ANALYSIS 

In  the  inviscid  program,  the  potential  flow  continuity  equation  takes 
the  form  of  the  Laplace  equation  in  terms  of  the  velocity  potential, 

V2^  =  0.  The  mathematical  problem  to  be  solved  requires  finding  a 
velocity  potential  which  is  harmonic,  and  satisfies  the  boundary  condi¬ 
tion  of  uniform  flow  at  infinity  and  flow  tangent  to  the  surface  of  the 


airfoil  and  the  centerline  of  the  jet.  Since  the  shape  of  the  jet  is 
initially  unknown,  Spence  (ref  10)  made  the  important  contribution  to 
solving  this  problem  by  relating  the  inertia  force  in  the  jet  to  the 
strength  of  an  equivalent  vortex  sheet, 


T/R  =  PU  y 


3-1 


In  this  equation  T  is  the  jet  thrust  (which  Spence  assumed  constant), 

R  is  the  local  radius  of  jet  curvature,  and  y  is  the  vortex  strength 
per  unit  length  along  the  jet. 

This  force  balance  becomes  an  additional,  dynamic  boundary  condition 
at  each  point  on  the  jet  sheet.  In  the  inviscid  analysis,  the  jet  thrust  is 
considered  to  vary  in  a  known  way.  This  variation  is  determined  from  the 
viscous  solution  for  the  turbulent  jet  computed  in  the  previous  itera¬ 
tion.  The  inviscid  solution  is  obtained  by  iterating  between  the  jet 
shape  and  the  jet  vorticity.  For  each  iteration,  the  problem  reduces 
to  solving  a  set  of  linear  influence  equations  with  unknown  vortex 
strengths.  The  boundary  condition  on  both  the  airfoil  and  jet  is  that 
of  no  net  normal  velocity  at  each  panel  control  point.  A  Gaussian 
elimination  solution  of  the  simultaneous  equation  is  used. 

The  airfoil  is  described  by  an  input  set  of  X  and  Y  coordinate  points. 
These  points  may  be  used  as  the  vortex  panel  endpoints,  or  an  increased 
number  of  panels  may  be  calculated  using  cosine  spacing  on  the  upper 
and  lower  surfaces  and  nonlinear  interpolation  of  the  input  points. 

An  airfoil  geometry  subroutine  performs  the  interpolations  if  required 
and  computes  the  individual  panel  angles  and  control  point  (panel 
centroid)  coordinates.  If  an  angle  of  attack  case  is  run,  all  geometry 
is  rotated,  and  all  displacements  are  included. 

In  the  solution,  the  airfoil  is  represented  by  distributions  of 
linearly  varying  vorticity.  The  paneling  scheme  begins  at  the  airfoil 
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trailing  edge  and  proceeds  forward  over  the  lower  surface  and  back  over 
the  upper  surface  to  the  trailing  edge.  Each  vortex  distribution,  with 
the  exception  of  the  first  consists  of  a  ramp-up  (increasing  from  zero) 
vorticity  over  one  airfoil  panel  and  a  ramp-down  (decreasing  to  zero) 
vorticity  over  the  next  airfoil  panel.  The  first  vortex  distribution 
uses  a  constant  vorticity  over  the  first  airfoil  panel  and  a  ramp  down 
over  the  second  airfoil  panel.  (This  distribution  was  used  to  better 
represent  the  near  constant  pressure  in  this  region  under  conditions  of 
large  jet  deflection  and  high  jet  momentum.)  This  paneling  scheme  gives 
a  number  of  unknown  vortex  strengths  which  is  one  less  than  the  number 
of  airfoil  panel  control  points,  at  which  the  boundary  conditions  were 
defined.  This  over-specified  problem  is  solved  by  the  method  of 
uniform  error,  (Bristow,  ref  19)  .  Figure  6  illustrates  the  airfoil 
vorticity  paneling  scheme.  The  jet  was  modeled  as  a  thin  sheet  of 
vorticity  using  the  assumptions  of  Spence  (ref  10).  Within  each  iteration, 
the  jet  shape  was  assumed  to  be  known.  For  the  first  iteration,  the 

Linear  Vorticity  Distribution 


N  Unknown  Vortex  Strengths, 

N  +  1  Control  Points 


Figure  6.  Vorticity  Distributions  on  the  Airfoil  and  Jet 
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empirical  jet  shape  equation  from  Chandra  (ref  20)  was  used.  For  sub¬ 
sequent  iterations,  the  jet  panel  angles  were  computed  using  an  itera¬ 
tion  scheme  described  in  a  later  section.  Briefly,  this  scheme  used 
the  previous  jet  shape  (panel  angles)  modified  as  guided  by  a  new  shape, 
computed  from  satisfying  the  dynamic  jet  boundary  condition,  Equation  3-1. 

Jet  panel  lengths  were  set  by  a  simple  geometric  progression  factor, 
usually  between  1.05  and  1.15.  Thus,  each  panel  length  was  larger  than 
the  one  just  upstream  by  the  factor  chosen.  The  first  jet  panel  length 
was  similarly  related  to  the  last  airfoil  upper  surface  panel  length. 

This  scheme  had  the  distinct  advantage  of  minimizing  numerical  problems 
associated  with  rapidly  changing  panel  sizes.  The  jet  panel  lengths 
and  angles  were  used  to  calculate  the  panel  endpoint  and  control  point 
(centroid)  coordinates  using  a  summation  beginning  at  the  airfoil 
trailing  edge. 

The  jet  was  represented  by  distributions  of  linearly  varying  vorticity. 
The  first  jet  vortex  was  a  single  ramp  down  function  over  the  first  panel. 
Subsequent  vortex  distributions  were  a  ramp  up  over  one  panel  and  a  ramp 
down  over  the  next  downstream  panel.  The  second  jet  vortex  distribution 
begins  on  the  first  jet  panel.  The  last  jet  vortex  distribution  ramps 
down  to  zero  strength  at  the  trailing  edge  of  the  last  jet  panel. 

When  required,  the  jet  entrainment  effects  are  represented  by  a 

number  of  sink  panels,  coinciding  with  the  jet  panels,  and  a  numerically 

-1/? 

integrated  sink  strength  which  decays  as  X  Ie-  from  the  jet  trailing  edge 
to  a  point  far  downstream  (typically  1000  chord  lengths).  The  sink 
strengths  are  computed  from  entrainment  output  from  the  viscous  program 
and  are  modeled  as  constant  over  each  jet  panel.  The  initial  value  of 
the  decaying  sink  function  is  set  equal  to  the  last  panel’s  constant 
strength.  The  effect  of  these  jet  sinks  is  to  modify  the  boundary  con¬ 
ditions  imposed  at  each  control  point.  The  jet  sink  panels  are 
illustrated  in  Figure  7. 
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Constant  Sink  Distribution 


Vortex  and  Sink  Influence  Equations 


The  equations  used  to  compute  the  velocities  induced  by  each  of  the 
linearly  varying  singularities  (vortices  or  sinks)  used  in  the  inviscid 
flow  analysis  are  given  below.  These  equations  were  developed  by 

Dvorak  and  Woodward  (ref  21).  Mote  that  each  equation  is  written  in  a 
coordinate  system  with  its  origin  at  the  left  hand  edge  of  the  singu¬ 
larity  panel  and  with  its  ordinate  axis  lying  along  the  singularity. 
Routines  within  the  inviscid  program  handle  rotation  and  translation  of 
an  influenced  point  into  the  panel  coordinate  system,  and  resolve  the 
induced  velocities  into  velocities  normal  and  tangential  to  the 
influenced  panel  control  point.  The  terms  in  the  equations  are  defined 
in  the  sketch  of  the  singularity  coordinate  system  sketched  in  Figure  8. 


The  velocities  induced  by  a  ramp  up  vortex  or  sink  distribution  are: 


U  w  JJl  log  \^-c)-l±  z_2  .  2i  rtan-l  1  .  tan'1  -S- 1  l 
v  s  2’Mc  s  -/  2  +  z2  c  [_  x  x-c  J  ; 
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Figure  8.  Singularity  Coordinate  System 


The  velocities  induced  by  a  constant  vortex  or  sink  distribution  are 


where  s  is  the  singularity  strength,  which  is  known  for  the  sink  dis 
tributions  and  set  equal  to  for  the  calculation  of  the  vortex 
influence  coefficients. 


The  velocities  induced  by  a  ramp  down  singularity  distribution  are 
calculated  by  subtracting  the  ramp  up  velocities  from  the  constant 
singularity  velocities.  Velocities  induced  by  the  composite  singu¬ 
larities  (ramp  up  on  one  panel  and  ramp  down  on  the  next  panel)  are 
simply  the  sum  of  the  two  contributions  computed  separately. 

Initial  Jet  Shape 

While  the  initial  jet  shape  is  not  critical  with  the  iteration  scheme 
used,  a  good  first  guess  minimized  computer  run  time.  For  the  present 
inviscid  flow  program,  the  initial  jet  shape  (in  terms  of  panel  angles) 
was  obtained  from  the  empirical  jet  angle  function  derived  by  Chandra 
(ref  20). 


4-  =  0.64  exp  (-.82  s/Cn)  +  0.36  exp  (-11.7  s/C^)  3-6 

ao 


For  any  initial  jet  deflection,  0o,  and  jet  momentum  coefficient, 

Op,  the  jet  angle  at  a  given  distance,  s,  along  the  jet  arc  could  be 
computed.  Jet  panel  lengths  were  computed  as  described  above,  and  the 
initial  panel  angles  were  then  the  average  of  the  angles  at  the  panel 
endpoints.  Initial  jet  position  was  computed  from  the  individual  panel 
angles  and  lengths. 

Solution  for  Singularity  Strengths 

The  airfoil  and  jet  panel  geometries  and  the  influence  equations,  3-2 
to  3-5  are  used  to  compute  the  influence  coefficient  matrix  and  the 
boundary  condition  vector.  Terms  in  the  influence  coefficient  matrix, 

Ai j ,  are  the  normal  velocity  induced  by  unit  strength  vortex  distribu¬ 
tion,  i,  at  the  control  point,  j.  The  terms  in  the  boundary  condition 
vector,  B j ,  are  equal  but  opposite  to  the  normal  velocities  which  exist 
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at  each  control  point  due  to  the  freestream  velocity,  V^,  and  the  jet 
sinks,  if  present.  The  following  matrix  equation  results: 


This  system  of  equations  is  solved  for  by  a  Gaussian  elimination 
procedure. 

Pressure  and  Force  Coefficients 


The  resulting  pressure  coefficients,  Cp^,  are  computed  from  the 
singularity  strengths  as  follows: 


On  the  airfoil: 


P  -  P 

r  local  „  2 

CP  =  ^ - 1  •  VT 
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where  V-p  is  the  total  tangential  velocity  at  the  control  point  due  to 
Va  and  all  singularities. 


On  the  jet: 


PU  -  pL 

CP  =  — r-  =  ’2Vt  r-vg 


3-0 


where  yavg  *s  t^le  average  of  the  vortex  strengths  at  the  panel  endpoints. 

Airfoil  and  jet  normal  force  and  axial  force  coefficients  are  computed 
by  summing  the  product  of  the  panel  pressure  coefficient  and  its 
respective  projected  area. 


Calculation  of  "New"  Pane!  Angles 

"New"  jet  panel  angle  refers  to  the  panel  angles  computed  from  the 
jet  vorticity  distribution  and  the  dynamic  boundary  condition, 

Equation  3-1.  They  are  different  from  the  next  iteration  jet  panel 
angles,  which  are  computed  from  the  last  iteration  angles  by  the 
iteration  scheme  described  below.  In  terms  of  the  local  jet  radius,  R, 
and  panel  length,  AS,  the  change  in  jet  flow  angle  from  one  end  of  a 
panel  to  the  other  is 

Afi  =  AS/R  3-10 

but,  from  Equation  3-1,  R  =  T/pVV.  Thus, 

A0  =  AS  pVy/T  3-11 

In  terms  of  the  local  pressure  coefficient,  Cp  =  />W/q  and  jet  thrust 
coefficient  Cp,  =  T/qc,  the  change  in  angle  is 

A0  =  -Cp  AS/CpC  3-12 

Each  "new"  jet  panel  angle  is  then  the  average  of  the  panel  endpoint 
angles  found  from  the  panel  pressure  coefficients,  local  momentum  co¬ 
efficient  (input  from  the  viscous  program)  and  panel  length. 

Iteration  Scheme 

An  iteration  scheme  was  written  which  eliminates  some  of  the  problems 
associated  with  ordinary  relaxation  techniques.  This  scheme  provides 
rapid  convergence  initially,  is  not  as  sensitive  to  oscillations  near 
the  jet  trailing  edge,  and  does  not  require  the  user  to  chose  a  weight¬ 
ing  factor. 
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This  iteration  scheme  does  not  use  the  "new"  jet  panel  angles 
directly.  Instead,  they  are  used  only  as  an  indication  of  the  direc¬ 
tion  of  change.  The  magnitude  of  the  change  in  panel  angle  is  deter¬ 
mined  by  the  factor  1  -  where  n  is  a  positive  integer  which  may  be 
different  for  each  panel.  Initially,  n  =  1  and  n  is  incremented  by  1 
each  time  the  required  direction  of  change  is  reversed.  If  the  "new" 
panel  angle  is  less  than  the  angle  used  in  the  previous  iteration,  the 
previous  angle  is  multiplied  by  the  factor  1  -  using  the  current 
value  of  n.  If  the  "new"  panel  angle  is  larger  than  the  angle  used 
in  the  previous  iteration,  the  previous  angle  is  divided  by  the  factor 
1  "  either  case,  the  next  iteration  jet  panel  angle  is  never 

allowed  to  change  by  an  amount  greater  than  "new"  angle  -  "old"  angle. 

The  iteration  continues  until  a  convergence  criterion  is  met.  The 
criterion  used  is  that  the  average  change  in  panel  angle  must  be  less 
than  0.01  times  the  initial  jet  deflection. 


TURBULENT  JET  ANALYSIS 


The  sink  strengths  which  represent  the  effect  of  the  jet  mixing  are 
calculated  from  a  viscous  solution  for  the  turbulent  jet.  It  is  pos¬ 
sible  to  compute  the  mixing  without  including  the  elliptic  effects  by 
taking  advantage  of  the  flow  geometry.  Since  there  is  a  primary  direc¬ 
tion  of  flow  (along  the  jet),  the  thin  shear  layer  approximation  can  be 
applied.  This  reduces  the  governing  elliptic  equations  to  a  parabolic 
set  which  can  be  solved  by  marching  along  the  jet  in  the  streamwise 
direction.  By  treating  the  radial  pressure  gradient  as  a  perturbation 
of  the  axial  pressure  gradient,  additional  simplifications  are  achieved. 
These  assumptions  result  in  considerable  savings  in  computer  storage  and 
running  time  compared  to  a  solution  of  the  complete  elliptic  problem. 
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The  governing  equations  are  derived  from  Reynolds  equations  for 
turbulent  flow  by  neglecting  streamwise  diffusion  and  upstream  convec¬ 
tion.  Thus,  the  effect  of  a  local  disturbance  is  transmitted  downstream 
by  pressure  and  convection  and  across  the  flow  by  pressure,  convection 
and  diffusion.  In  cylindrical  coordinates,  (r,  9),  the  equations  take  . 
the  general  forms 

Conservation  of  Mass 


au  av  v 

r  69  dr  r 


=  0 
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Streamwise  Momentum 


au  au  uv  __  i  ap_  _  auV  2u'v' 

r d9  dr  *  r  ~  ~  P  rd 9  ~  dr  "  r 


Radial  Momentum 
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av  v  av  _  u£ 

rd 9  dr  r 


.  I  J2.(p 
P  dr  ' 
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in  which  U  and  V  are  the  streamwise  and  radial  velocity  components.  The 
radial  perturbation  pressure  p'(r,  9)  is  decoupled  from  the  mean  stream- 
wise  pressure,  P (9). 

These  equations  do  not  describe  the  effect  of  curvature  on  the  flow. 
The  "extra"  terms  originate  in  the  choice  of  coordinate  system  and  do 
not  imply  any  additional  physical  processes.  The  Coriolis  acceleration, 
UV/r  =  «»V  and  the  centrifugal  acceleration  uVr  =  «U,  arise  from  rotation 
of  the  velocity  vector  with  respect  to  the  coordinate  system.  Similarly, 
the  extra  stress  term  is  a  consequence  of  having  chosen  a  cylindrical 
volume  element;  that  is,  since  the  sides  of  the  element  are  not  parallel, 
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the  shear  stress  has  a  component,  u'v’/r,  in  the  streamwise  direction. 
Thus,  the  use  of  cylindrical  coordinates  is  a  convenience  in  this  case, 
and  the  additional  terms  merely  account  for  the  difference  between  the 
actual  accelerations  of  the  fluid  and  those  which  appear  to  take  place 
in  the  curvilinear  coordinate  system. 


To  facilitate  solution,  the  governing  equations  are  transformed  into 
anu-x  coordinate  system  as  shown  in  Figure  9.  This  provides  the 
advantage  of  expanding  the  grid  point  spacing  as  the  jet  spreads.  Thus, 
the  maximum  number  of  grid  points  are  always  contained  within  the  jet, 
and  there  is  no  need  to  include  additional,  unnecessary  grid  points  for 
the  external  stream.  The  stream  function,  is  defined  as 
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Figure  9.  Transformed  Coordinate  System  for  Turbulent  Jet 
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The  variable,  <#,  is  a  dimensionless  stream  function  defined  by: 


*E  -  *i 
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The  ^  coordinate  lines  are  not  true  stream  lines.  They  are  used  as  a 
convenient  method  to  handle  flow  in  the  cross  stream  direction. 


Turbulence  Equations 


The  effect  of  curvature  is  to  increase  the  turbulent  mixing  rate. 
Recent  proposals  to  model  this  effect  have  been  based  on  the  turbulence 
kinetic  energy  equation.  According  to  the  usual  eddy  viscosity  assump¬ 
tion,  the  turbulent  stress  is  first  expressed  in  terms  of  a  turbulent 
viscosity,  Ht,  and  the  mean  strain  rates 


u’v'  =  4,. 


,6U  0. 

or  ’  r 
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Again,  no  new  physical  process  is  implied  by  the  extra  strain  rate, 

U/r.  Following  Launder  and  Spalding  (ref  22),  the  turbulent  viscosity 
is  assumed  to  depend  on  two  parameters:  the  turbulence  kinetic  energy, 
k,  and  its  rate  of  dissipation,  €.  From  dimensional  arguments,  the 
expression  for  the  turbulent  viscosity  is 


^t 


.2 

Vk 

€ 
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in  which  c^  is  an  empirical  constant.  A  procedure  for  modifying  the 
turbulence  kinetic  energy  equations  to  account  for  the  effect  of  curva¬ 
ture  has  been  developed  by  Schum,  et  al.  (ref  23).  This  method  has  been 
used  here.  Briefly,  the  effect  of  curvature  is  to  transfer  energy  from 
the  streamwise  component  of  the  turbulence  to  the  transverse  component. 


Since  it  is  the  transverse  component  which  produces  most  of  the  mixing 
in  thin  shear  layers  (Bradshaw,  ref  24),  this  transfer  of  energy  increases 
the  mixing. 

However,  curvature  has  no  effect  on  the  total  turbulence  energy, 
because  energy  is  simply  transferred  from  one  component  to  another. 

Thus,  the  turbulence  kinetic  energy  equation  cannot  be  used  to  calculate 
the  effects  of  curvature.  An  additional  production  term  was  therefore 
added  to  the  total  turbulence  energy  equation  to  represent  the  increase 
of  the  transverse  energy  component  by  the  curvature.  In  cylindrical 
coordinates,  the  equation  for  the  turbulence  energy  becomes 


e&- 

dt 


,dU  U 


<§*  -  *)  -  pf  +  -2-  (r* 


dr 


dr  <*k  dr 


2U 

r 
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in  which  the  first  term  on  the  right  hand  side  represents  the  usual 
production  mechanism,  the  second  term  represents  dissipation,  the 
third  term  represents  diffusion,  and  the  last  term  represents  the 
effect  of  curvature. 


An  equation  derived  by  Chambers  and  Wilcox  (ref  25)  was  used  to  compute 
€,  the  dissipation  of  the  turbulence  energy.  It  contains  the  same 
curvature  correction  as  the  turbulence  energy  equation: 


P 


d€ 

dt 


(C1g  -  c2p£)  k  +  ^(^rlf) 


+  U  '  v '  — 


2U 

r 
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This  turbulence  model  has  five  empirical  constants:  cp,  cj.,  C2,  <t^» 
and  <*£.  The  values  of  these  constants  recommended  by  Launder  and 
Spalding  (ref  22)  were  used  for  this  analysis. 
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Finite  Difference  Equations 


The  governing  mean  flow  and  turbulence  equations  are  solved  simul¬ 
taneously  by  iteration.  In  order  to  simplify  the  program  logic,  the 
momentum  and  turbulence  equations  are  expressed  in  the  general  form: 


,  .  .  s  d  d±s 

(a  +  b->  eZ  =a^r  (rgT> 


in  which 


mi 

3  ‘  *E  - 


and 
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The  general  dependent  variable  ^  represents  either  U,  V,  k,  or  e, 
while  T  and  S  represent  the  corresponding  diffusion  and  source  terms. 
Also,  m  represents  the  mass  flux  across  the  jet  boundaries  shown  in 
Figure  9.  With  all  the  equations  expressed  in  this  same  form,  the 
same  solution  procedure  can  be  used. 

A  marching  solution  is  performed  in  which  the  task  is  to  obtain 
values  of  the  flow  variables  at  a  downstream  station  when  the  values 
at  the  upstream  station  are  known.  With  a  given  set  of  initial  con¬ 
ditions,  one  can  then  make  a  series  of  forward  steps  in  the  streamwise 
direction  to  obtain  the  solution  for  the  entire  flow  field.  The  flow 
field  is  modeled  by  specifying  a  number  of  control  volumes  across  the 
flow  in  a  direction  normal  to  the  mainstream  direction,  x,  with  grid 
points  centered  at  the  control  volume  faces  upstream  and  downstream 
(Figure  10). 

The  method  of  specifying  control  volumes  in  the  program  is  one  where 
the  control  volume  boundaries  are  input  in  the  «  coordinate  system  (at 
x  ■  0)  and  the  respective  grid  points  are  then  located  in  the  center  of 
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Figure  10,  Control  Volume  and  Coordinate  System 

each  control  volume  (midway  between  y+  and  y~).  The  points  Y+  and  Y" 
represent  adjacent  grid  points,  above  and  below  grid  point  P<j.  The 
values  of  Ay  can  change  in  the  x  direction,  depending  on  the  values  of 
the  U  velocities  in  the  control  volume.  The  finite-difference  equations 
are  then  formed  by  integrating  the  governing  partial  differential  equa¬ 
tions  over  the  control  volumes.  In  order  to  express  the  various  terms 
in  the  equations  in  algebraic  form,  the  following  practices  are  used: 

1)  The  value  of  the  variable  at  the  central  grid  point  is  assumed 
to  prevail  over  the  entire  face  of  the  control  volume. 

2)  The  downstream  values  are  considered  to  prevail  over  the  entire 
forward  step. 

3)  The  dependent  variables  are  assumed  to  vary  linearly  between  the 
grid  points  (this  is  modified  when  convection  is  large). 
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4)  Strict  compatibility  is  maintained  between  the  neighboring 

control  volumes,  ie.,  the  flux  leaving  a  control  volume  through 
a  face  is  made  exactly  equal  to  the  flux  entering  the  adjacent 
control  volume  through  the  same  face. 

The  resulting  algebraic  equations  (finite-difference  equations)  for 
each  grid  point  have  the  following  general  form: 

Ap  <^p(j  =  Ayt-  ^y+  Ay-  ^y-  +  Apu  ^Pu  ®  3-24 

where  Ay+  and  Ay-  represent  the  effect  of  cross  stream  convection  and 
diffusion;  Apu  represents  the  upstream  convection;  B  is  the  constant 
part  of  the  source  term 


S  =  B  +  Sp  <t> 
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and  Ap  is  defined  by 


Ap  =  Ay+  +  Ay-  +  Apu  -  Sp  3-26 

The  symbolic  term,  is  used  to  represent  the  U,  k,  or  €  parameters 
at  grid  locations,  Y*,  Y" ,  u  or  d.  Table  I  presents  the  corresponding 
source  terms.  The  value  of  the  A  coefficients  depend  on  which  parameter 
is  being  considered.  The  coefficients  Ay+  and  Ay-  in  Equation  3-24  are 
defined  as: 


aY+  =  d£+  *  2  S+ 
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Ay-  =  D*.  +  |  Cy-  3-28 

D  represents  diffusion  and  C  represents  convection  through  the  y+  or  y~ 
boundary.  The  subscripts  for  D  and  C  indicate  that  these  quanti  are 
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where  r  ,  <«,  and  Ax  are  the  values  of  the  exchange  coefficient, 
dimensionless  stream  function,  and  x  stepsize,  respectively,  of  the 
control  volume  boundary  being  considered.  Ay'  is  the  distance  between 
the  two  grid  points  on  either  side  of  the  control  volume  boundary. 

The  coefficient  Apu  represents  the  convection  from  upstream  and  is 
calculated  by 


APu  =  -  *I>U  A“  3-31 

where  tota^-  upstream  flow  rate  and  Au>  is  the  change 

in  <a  across  the  control  volume  being  considered. 


Once  the  solution  of  the  differential  equations  is  complete  and 
values  of  the  U  velocities  and  AP  have  been  calculated,  the  locations 
of  the  grid  points  and  control  volume  boundaries  in  terms  of  y  must  be 
found.  This  is  done  by  using  the  definitions  of  and  »  to  obtain 


dy 


(»E  - 
PU 


d  u 
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Now  by  using  this  expression  and  the  known  u>  values  that  are  input  for 
the  control  volume  boundaries  the  corresponding  Ay  can  be  found  from 


(K  -  *I> 

■  -fer- 


(An). 


n  =  2 ,3 , . . .  ,N 
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for  each  n  control  volume.  The  y  values  of  the  control  volume 
boundaries,  y+,  are  then  easily  iound  by  knowing  the  boundary  value, 
yj»  (which  can  be  set  to  0)  and  the  (Ay)  values: 


<»E  -  ^  + 

(pU)n  y  n-1 


n  =  2,3,...,  N+l  3-34 
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will  equal  y£. 


where  y+i  =  Yi  and  y+N+1 

The  y  location  of  the  grid  point  n  is  calculated  by  using  the 
centered  grid  point  approach  so 

yn  =  \  (y+n  +  y4^)  n  =  2,3,...,  N=1  3-35 

Pressure  Equations  for  Curved  Flow 

In  the  viscous  program,  two  pressure  gradients  are  considered--one 
in  the  streamwise  direction  and  the  other  in  the  cross  stream  direction, 
a  result  of  jet  curvature.  With  the  addition  of  curvature  effects,  the 
change  in  the  local  pressure  due  to  the  centrifugal  forces  must  be 
included  in  the  mainstream  momentum  equation.  However,  in  parabolic 
flow,  these  pressure  changes  must  not  change  the  total  integrated  mean 
pressure  difference  times  area  term  for  a  cross  section.  In  other 
words,  only  one  value  of  mean  pressure  exists  at  a  particular  x  station. 
The  local  pressure  variation  due  to  curvature  is  calculated  by  creating 
an  additional  pressure  value,  P',  from 


dP'  PU2 

d  y  r 
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where  r  is  the  local  curvature.  Upstream  values  of  the  velocity,  Uu, 
are  used  for  downstream  values  since  the  latter  have  yet  to  be  deter¬ 
mined.  With  small  Ax  steps,  as  used  here,  this  is  a  good  approximation. 


For  the  mean  or  global  pressure  to  have  only  one  value  at  a  particular 
x  station,  it  is  necessary  that: 

I  Pj  dAj  =  0  j  =  1,  2,...  N  3-37 
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where  N  is  the  number  of  control  volumes.  The  values  of  Pj  are  used 
to  determine  a  gradient  dP'/dx  for  use  in  the  source  term  in  the  U 
momentum  equation  (Table  I). 

Boundary  Conditions 

The  viscous  program  requires  boundary  conditions  to  be  specified  at 
the  inlet  plane  and  along  the  jet  surfaces.  At  the  inlet  plane  (nozzle 
exit)  a  top  hat  velocity  profile  is  used.  The  radial  perturbation 
pressures  are  determined  from  Equation  3-36  in  which  the  radius  of 
curvature  is  specified  by  the  inviscid  program.  The  U  velocity  at  the 
inlet  is  derived  from  a  stream  tube  analysis,  Figure  11.  By  combining 
the  Bernoulli  and  continuity  equations  with  the  definition  of  C^,  the 
initial  jet  velocity,  Uj^  can  be  obtained. 


Figure  11.  Stream  Tube  for  the  Deflected  Jet 
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From  Bernoulli 


Pi  +  2  pUJ2,i  =  P®  +  2pUJ2,, 


and  continuity 


P  U j  ,  A-  =  P  Uj  A» 

J  ,  1  1  J  ,00 


It  can  be  shown  that  U.  .  can  be  obtained  from  the  quadratic 

J ,  l 


2  2(Pj.  -  ^)U  2  /U«2  C  C\2 


At  any  downstream  distance,  x,  the  velocity  profile  is  dependent  on 
the  solution  of  the  momentum  equations  with  prescribed  velocities  as 
boundary  conditions  at  the  inner  and  outer  surface  of  the  jet.  There 
are  two  ways  to  obtain  the  edge  velocities,  one  is  to  use  the  velocities 
from  the  inviscid  program  and  the  other  is  to  use  the  upstream  radial  pressures 
from  the  viscous  program.  The  latter  approach  was  used  in  the  iteration 
procedure  involving  the  inviscid  and  viscous  programs.  This  serves  as 
a  check  on  the  degree  of  convergence  between  the  two  programs  since  the 
velocities  at  corresponding  locations  should  be  the  same  when  convergence 
is  achieved.  The  inner  and  outer  jet  surface  velocities  are  determined 


from: 


<Ul)Viscous  =  VI  <P  +  Pi> 


^Viscous  "I  7  <*  +  **> 


in  which  P'  and  P2.  are  determined  from  the  previous  upstream  radial  pres- 
I  b 

sure  distribution  (Eq.  3-36),  calculated  in  the  viscous  program.  The  mean 
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pressure,  P,  is  obtained  from  the  inviscid  velocities. 
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The  velocity  profiles  in  the  jet  are  also  dependent  upon  the  entrain¬ 
ment  at  the  inner  and  outer  jet  surface,  mg  and  thj.  The  values  are 
obtained  from  the  relation  given  by  Patankar  and  Spalding  (ref  18) 


■©). 


jet  edge 


in  which  the  viscosity  is  calculated  in  the  solution  of  the  momentum  and 
turbulence  equations.  In  some  cases,  step-to-step  oscillations  in  the 
width  of  the  jet  can  occur.  One  of  the  many  approaches  to  minimize  these 
oscillations  (Patankar  and  Spalding,  raf  18)  is  to  utilize  a  multiplication 
factor,  F,  in  the  entrainment  calculations,  such  that 
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where  F  is  altered  from  one  Ax  step  to  the  next  by 
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where  jUfH-l  “  Ufj)  *-s  the  absolute  value  of  the  difference  in  velocity 
between  the  boundary  and  adjacent  grid  point,  and  AU*  is  the  velocity 
difference  specified  by  the  user  as  the  minimum  to  which  he  wants  the 
difference  between  Ufj+i  and  1%  to  fall.  The  exponent  controls  the  rate 
of  change  of  thickness  if  oscillations  are  a  problem.  Setting  it  equal 
to  0.05  has  generally  worked  well.  AU*  is  generally  specified  as  a 


percentage  of  the  maximum  change  in  velocity  in  the  velocity  profile, 
say  one  or  two  percent  of  Umax  “  umin  •  Even  though  the  value  of  F 
may  vary  from  step  to  step,  on  the  average,  the  correct  entrainment  is 
calculated. 

Solution  of  Finite  Difference  Equations 

For  each  grid  node  (excluding  boundary  nodes)  there  is  a  correspond¬ 
ing  finite  difference  equation  (Equation  3-24)  representing  the  momentum 
and  turbulence  equations.  If  there  are  M  nodes  used,  there  are  (M  -  2) 
equations  to  be  solved.  These  finite-difference  equations  are  solved 
by  using  the  tri-diagonal  matrix  algorithm  having  the  generalized  format 

Aj  *j  =  Bj  *j+l  +  +  D*  3-47 

for  j  =  2,3,...,M-1.  <t>i  and  are  prescribed  boundary  conditions,  such 
as  the  jet  edge  velocities.  Equation  3-47  represents  Equation  3-24  in 
which  the  last  two  terms  have  been  combined  into  the  single  value,  D*. 

When  calculating  the  coefficients  for  the  equations,  several  values 
of  flow  variables  and  distances  are  needed.  Upstream  values  are  used 
when  downstream  values  are  not  available.  No  iterations  per  forward 
step  are  performed.  These  iterations  could  improve  the  solution;  how¬ 
ever,  the  improvement  is  insufficient  to  warrant  the  increased  work  for 
two-dimensional  parabolic  problems.  This  lack  of  iteration  is  only  a 
problem  if  the  forward  steps  become  too  large  and  one  is  covering  an 
area  where  the  flow  variables  change  a  lot  over  one  step.  For  the  cases 
studied  here,  the  forward  Ax  step  was  only  0.0005  ft  over  a  total  jet 
length  of  about  10  ft.  Changes  in  the  flow  variables  over  such  a  small 
distance  as  0.0005  ft  are  relatively  small. 
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The  only  readily  available  data  for  a  turbulent  jet  in  a  coflowing 
stream  with  a  streamwise  pressure  gradient  is  that  of  Fekete  (ref  26). 
Tests  were  conducted  so  as  to  provide  self  preservation.  According  to 
Newman  (ref  27),  self  preservation  is  possible  only  if: 


Ho 

U* 


constant 
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and  U*  a  X  m 

where 

+  Uq/u* 

m  =  -  i — 

3  V  2  +  2  Uq/Uod 
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Measured  velocity  profiles  of  Fekete  showed  that  self  preservation  was 
achieved.  Velocity  profiles  were  also  calculated  with  the  viscous 
program  for  a  set  of  conditions  identical  to  Fekete,  namely  for 
uo/ua>  ■  0.95.  Analytically  predicted  velocity  profiles  also  had 
similarity  and  are  shown  in  Figure  12  to  be  in  good  agreement  with 
measured  values. 


Experimental  values  of  the  excess  momentum  were  obtained  by  integrat¬ 
ing  Fekete's  similarity  solution  for  the  jet  velocity  profiles. 

Computed  values  of  the  excess  momentum  were  within  six  percent  of  the 
measured  values.  The  computed  jet  spreading  rate  of  0.037  agreed 
favorably  with  the  measured  value  of  0.036,  and  the  computed  gradient 
of  the  excess  momentum  was  within  27.  of  the  measured  values.  Since  the 
turbulence  constants  were  not  adjusted,  this  is  very  good  agreement. 
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Figure  12.  Comparison  of  Calculated  and  Measured  Jet 
Ve locities 

MATCHING  PROCEDURE 

The  analysis  of  the  jet  flap  is  matched  to  the  analysis  of  the 
turbulent  jet  by  iterating  between  solutions  for  the  shape  of  the  jet 
flap  and  the  variation  of  the  jet  thrust,  as  follows:  for  given  dis¬ 
tributions  of  jet  thrust  and  sink  strength,  the  inviscid  jet  flap 
program  yields  an  estimate  of  the  jet  shape.  The  curvature  and  pres¬ 
sure  distribution  along  the  jet  sheet  then  provides  a  set  of  boundary 
conditions  for  calculating  the  jet  thrust  and  sink  strengths  with  the 
turbulent  mixing  program.  These  solutions  are  iterated  until  the 
computed  reduction  of  jet  thrust  and  the  jet  drag  induced  on  the 
airfoil  converge  to  within  acceptable  limits.  The  iterations  are 
started  with  a  classical  jet  flap  solution  in  which  the  thrust 
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of  the  jet  is  assumed  to  be  constant,  and  the  sink  strengths  are  all 
zero.  Generally,  three  to  four  iterations  are  required  to  obtain  con¬ 
vergence. 


The  thrust  and  momentum  of  the  jet  are  evaluated  in  the  following 
way:  The  net  thrust  on  the  airfoil  can  be  calculated  by  integrating 
the  computed  velocity  profile  at  any  point  in  the  wake.  For  simplicity, 
consider  the  case  in  which  the  jet  is  undeflected,  as  shown  in  Figure 
13.  If  the  upstream  control  surface  is  far  enough  away  for  the  flow 
across  it  to  be  undisturbed,  the  momentum  equation  yields: 


T  = 


(P  +  pU2)dy  - 


(Poo  +  PU2)dy 


3-51 


By  factoring  and  recombining  the  terms,  we  get 

+00  +00  +*» 


T  -  p 


/ 


UOJ-UJdy  +  pU* 


/ 


(U-Uco)dy  + 


/ 


(P-PJdy 
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The  continuity  equation  expresses  the  fact  that  the  difference  in  mass 
flow  between  the  upstream  and  downstream  stations  is  equal  to  the  mass 
added  by  the  jet;  that  is. 


T 

/ 


(U  -  Uoo)  dy 
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Substituting  this  relation  into  the  momentum  equation  then  yields  the 
following  expression  for  the  net  thrust  on  the  airfoil: 


T 


®  j  Uoo  +  P 


+>» 


U(U  -  UcJdy  + 


(P  -  P00)dy 


3-54 
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Figure  13.  Control  Surfaces  for  Evaluating  Jet  Thrust 

The  first  term  is  the  source  thrust  generated  by  the  added  mass  of  the 
jet.  For  turbojet  engines,  which  draw  the  initial  mass  of  the  jet  from 
the  atmosphere  (except  for  a  very  small  fuel  fraction)  the  source  thrust 
is  negligably  small.  This  term  is  also  zero  in  a  pure  wake.  On  the 
other  hand,  for  rocket  engines,  which  generate  the  initial  jet  mass  from 
internal  fuel,  this  term  is  significant.  It  must  also  be  included  in 
evaluating  the  thrust  of  the  jet  flap  model  since  the  jet  mass  is  added 
to  the  tunnel  mass  flow  from  inside  the  model.  The  second  term  is  the 
integral  of  the  excess  momentum;  it  represents  the  part  of  the  initial 
jet  momentum  which  is  greater  than  the  displaced  free  stream  momentum. 

In  a  pure  wake,  it  would  represent  the  drag  of  the  body.  The  last  term 
represents  the  pressure  force  due  to  the  perturbations  of  the  velocity 
field.  Close  to  the  model,  this  pressure  term  is  significant.  However, 
it  vanishes  at  points  far  downstream.  As  it  does,  the  excess  momentum 
increases,  because  the  net  thrust  is  the  same  at  every  station  behind 
the  model. 


is  used  in  the  inviscid  program  for  calculating  the  jet  shape.  Us  is 
the  velocity  in  the  main  stream  on  either  side  of  the  jet.  An  average 
of  the  values  at  the  jet  boundaries  is  used  for  the  value  of  Us  within 
the  jet.  Outside  the  jet,  where  U  =  Us,  the  integral  vanishes.  At  down- 
stram  infinity,  where  Us  =  U^,  the  jet  momentum  flux  does  equal  the  net 
thrust  of  the  jet. 


The  sink  distribution  represents  the  drag  induced  on  the  airfoil  by 
the  jet.  In  principle,  the  sinks  can  be  determined  from  the  streamwise 
variation  of  the  excess  mass  within  the  jet  (The  variation  of  the  excess 
mass  must  be  used,  rather  than  the  total  entrainment,  because  only  the 
change  in  excess  mass  affects  the  flow  outside  the  jet.)  However,  due 
to  the  arbitrary  way  in  which  the  jet  boundary  is  defined,  it  is  more 
accurate  to  evaluate  the  sink  strengths  directly  from  the  change  in  jet 
thrust. 

The  force  on  a  section  of  the  jet  and  the  equal  but  opposite  reaction 
on  the  airfoil  are  shown  in  Figure  14.  The  reaction  force  on  the  airfoil 
has  two  components:  pUqr ,  a  component  due  to  the  interaction  between  the 
sink  and  airfoil  singularities,  and  pU^dr,  a  component  due  to  the  change 
which  the  sink  produces  in  the  strength  of  the  airfoil  singularities. 
These  are  also  shown  in  the  figure.  Since  the  horizontal  force  induced 
by  the  sink  must  equal  the  horizontal  component  of  the  thrust  loss,  an 
estimate  of  the  sink  strength  at  each  point  on  the  jet  is  obtained  from 
the  relation 
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Figure  14.  The  Equal  but  Opposite  Drag  Forces  on  the  Airfoil  and  Jet 

AT*  cos  8  =  pUqT  cos  X  3-56 

in  which  AT*  is  the  local  change  in  the  actual  thrust  of  the  jet,  and 
pUgFcosA  is  the  net  horizontal  force  induced  on  all  the  other  singularities 
by  the  sink.  Alternately,  the  equal  but  opposite  force  induced  on  the 
sink  by  all  the  other  singularities  may  be  used  on  the  right  hand  side. 

In  the  next  inviscid  calculation,  the  strength  of  the  airfoil  singularities 
are  corrected  to  include  the  influence  of  the  sink  distribution.  There¬ 
fore,  when  the  iterations  have  converged,  both  components  of  the  force 
on  the  airfoil  are  determined. 

The  local  thrust  of  the  jet,  T*,  is  defined  to  be  the  thrust  obtained 
by  an  isentropic  expansion  (without  additional  mixing)  of  the  local  mass 
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of  the  jet  to  the  pressure  at  infinity.  It  is  necessary  to  use  the  local 
thrust  for  defining  the  sink  strengths,  rather  than  the  local  excess 
momentum  because  the  excess  momentum  does  not  include  the  contribution  of 
the  pressure  term  to  the  jet  thrust. 

The  local  thrust  of  the  jet  is  evaluated  from  the  local  jet  velocities 
in  the  following  way.  The  local  thrust  is  defined  as 

T*  =  mjU^  +  p  J"  U*(U*  -  U«Jdy  3-57 


in  which  U*  prescribes  the  velocity  distribution  of  the  expanded  jet. 
Since  the  expansion  is  assumed  to  occur  without  mixing,  the  continuity 
equation  for  each  stream  tube  yields 


pU*dy*  =  pUdy 


3-5fi 


Thus,  the  local  thrust  can  be  written 


T*.  =  mjU  +  p 


+=> 

/‘ 


U(U*  -  Ujdy 
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Then,  since  U*  is  defined  by  Bernoulli's  equation, 

j  PU*2  +  P*,  =  j  pU2  +  P  3-60 

the  local  thrust  can  be  expressed  in  terms  of  the  local  jet  velocity 
distribution: 


T*  =  +  p  j  U([U2  +  2  (P  -  Pj/p]15  -  lUdy 


3-61 


In  the  irrotational  flow  on  each  side  of  the  jet,  U  =  Us  and 


2 (P  -  P^/p  =  U2  -  U2 
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(Us  has  a  different  value  on  each  side  of  the  jet.)  As  a  result,  this 
integral  differs  from  zero  only  within  the  jet  so  that  the  local  thrust 
can  be  evaluated  from  the  computed  velocity  profiles  within  the  jet. 


As  a  footnote  to  this  discussion,  Jones  (Schlichting,  ref  23)  derived 
Equation  3-61  for  use  in  determining  airfoil  drag  from  a  wake  survey. 
However,  by  expanding  the  radical  in  powers  of  2A P/pU2,  it  can  be  seen 
that  this  equation  does  not  give  the  drag  (or  thrust)  on  the  airfoil: 


This  series  converges  for  AP  <  pU2,  but  is  equal  to  the  actual  drag  only 
when  AP  <<  pU2.  Thus,  Jones'  method  does  not  yield  the  actual  drag  when 
AP  -  pU2.  His  error  was  the  assumption  that  evaluating  U*  with  Bernoulli’s 
equation  gives  the  actual  drag.  To  obtain  the  drag,  U*  must  be  defined 
from  the  drag  equation 


which  yields 


U*  =  U(1  +  AP/pU2) 


3-65 


On  the  other  hand,  the  integrand  does  not  go  to  zero  outside  the  jet  with 
this  definition  of  U*,  so  that  there  is  no  advantage  over  using  Equation 
3-54.  We  have  therefore  used  the  method  of  Betz  (Schlichting,  ref  2S)  to 
measure  the  jet  thrust  in  the  experimental  phase  of  this  study. 


Jet  Drag  of  a  Cylinder 

In  order  to  illustrate  the  results  of  this  matching  technique,  the 
jet  drag  induced  on  a  cylinder  has  been  calculated.  A  total  of  36 
vortex  panels  were  used  to  represent  the  cylinder,  and  100  sink  panels 
were  used  for  the  jet.  The  jet  was  considered  to  originate  at  the  rear 
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stagnation  point  and  trail  straight  back  along  the  stagnation  streamline 
A  Reynolds  number  of  Re  =  500,000  and  jet  thrust  coefficient  of  Cu  =  1 
were  used.  However,  it  was  assumed  that  the  viscous  drag  of  the 
cylinder  was  zero  and  that  there  was  no  separated  wake. 

A  total  of  three  iterations  were  required  for  convergence.  In 
Figure  15  the  computed  surface  pressure  distribution  for  this  case  is 
compared  to  the  computed  distribution  for  the  classical  case  without  a 
jet.  It  can  be  seen  that  the  effect  of  the  jet  is  to  accelerate  the 
flow  around  the  cylinder  and  reduce  the  base  pressure.  This  produces 
a  net  pressure  drag  of  magnitude  D/Tj  =  .207  on  the  cylinder. 

The  streamwise  variation  of  the  excess  momentum  in  the  jet  is  shown 
in  Figure  16.  At  downstream  infinity,  the  excess  momentum  plus  the 
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Figure  15.  Comparison  of  Cylinder  Surface  Pressure  Distributions 
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Figure  16.  Variation  of  Excess  Momentum  in  the  Jet 

source  thrust  give  the  net  force  on  the  cylinder.  By  X/D  =  7,  the 
reduction  in  jet  thrust  is  AT/Tj  =  .189.  Considering  the  difficulty 
in  accurately  determining  the  drag  from  surface  pressure  integrations, 
this  was  felt  to  be  good  agreement,  and  additional  iterations  were  not 
performed.  A  solution  was  also  attempted  using  the  sink  strengths 
given  by  the  total  jet  entrainment,  but  the  jet  drag  on  the  cylinder 
approached  a  value  near  D/Tj  =  0.5,  which  was  roughly  twice  the  change 
in  jet  thrust. 
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IV.  EXPERIMENTAL  METHODS 


WIND  TUNNEL  FACILITY 

Tunnel  Description 

The  wind  tunnel  in  which  the  jet  flap  model  was  tested  is  a  two- 
dimensional  tunnel  designed  especially  for  powered  lift  testing.  A 
photograph  of  this  tunnel  is  shown  in  Figure  17.  The  tunnel  is  an 
open  circuit  return  type  driven  by  a  downstream  multi-tube  ejector. 

The  test  section  is  51  cm  wide  and  366  cm  high.  With  the  20  cm 
chord  model  of  this  test,  the  height-to-chord  ratio  for  this  tunnel 
is  18.0.  This  ratio  minimizes  wall  effects  and  allows  sufficient 
space  for  the  jet  to  develop  downstream  for  the  large  jet  deflections 
and  momentum  coefficients  tested.  For  example,  with  Cy  -  50°,  the 
velocity  overspeed  is  estimated  to  be  1.001  U®-  The  test  section 
extends  approximately  8  chord  lengths  (158  cm)  upstream  and  12  chord 
lengths  (23  cm)  downstream  of  the  model  quarter  chord.  A  smooth 
entrance  to  the  test  section  is  provided  by  a  61  cm  long  elliptical 
bellmouth.  The  vertical  walls  of  the  tunnel  diverge  from  51  cm  apart 
at  the  downstream  end  of  the  bellmouth  to  55  cm  apart  at  the  tunnel 
exit  to  allow  for  wall  boundary  layer  growth.  The  model  was  mounted 
in  the  test  section  213  cm  above  the  tunnel  floor.  The  tunnel  was 
constructed  of  plywood  and  lumber. 

A  multi-tube  ejector  was  located  at  the  downstream  end  of  the  test 
section  to  power  the  tunnel.  The  ejector  primary  nozzle  arrangement 
consisted  of  304  separate  .457  cm  I.D.  nozzles  uniformly  spaced  in  the 
ejector  plane.  These  nozzles  were  fed  by  two  10  cm  diameter  steel 
plenums,  located  outside  the  tunnel,  one  on  either  side.  The  ejector 
gives  an  entrained  to  primary  mass  flow  ratio  of  20.  With  the  air 
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supply  available  for  this  test,  the  maximum  tunnel  speed  was  approxi¬ 
mately  24  meters  per  second. 

The  model  was  mounted  in  aluminum  plates  set  flush  in  the  tunnel's 
wooden  walls.  These  plates  contained  circular  insets  which  allow  the 
model  to  be  rotated  to  angles  of  attack,  between  -10  and  +20  degrees. 
Windows  were  provided  in  the  inserts  to  review  the  model  during  testing. 
A  photograph  of  the  aluminum  plates  is  shown  in  Figure  18. 


Figure  18.  Airfoil  Model  Endplates 


Two  boundary  layer  blowers  were  mounted  on  each  test  section  side- 
wall.  A  30  cm  long  BLC  blower  was  located  7.6  cm  upstream  of  the  model 
leading  edge  and  was  set  at  an  angle  of  15  degrees  from  the  vertical. 

A  15  cm  long  blower  was  attached  to  the  model  mount  inserts,  with  the 
nozzle  located  at  the  midchord  blowing  parallel  to  the  chordline.  All 
four  blowers  were  fed  from  a  common  supply  line  and  operated  at  the 
same  supply  pressure.  The  blower  nozzle  gaps  were  set  to  .100  cm. 

Tunnel  Calibration 

Calibration  of  the  wind  tunnel  test  section  was  performed  in  two 
steps.  The  first  stage  was  to  adjust  the  vertical  wall  divergence  to 
obtain  a  constant  static  pressure  distribution  along  the  entire  length 
of  the  test  section.  Then  the  total  and  static  pressure  distributions 
in  the  area  where  the  model  would  be  located  were  surveyed.  This  survey 
data  was  used  to  generate  a  calibration  curve  of  test  section  conditions 
versus  a  measured  sidewall  surface  static  pressure.  The  results  of  the 
calibration  procedures  are  discussed  below. 

A  5  cm  diameter  pipe  containing  static  pressure  taps  every  15.25  cm 
was  mounted  equidistant  between  the  tunnel  walls  and  213  cm  above  the 
tunnel  floor  (at  the  future  model  chord  line) .  The  first  static  tap 
was  at  the  downstream  end  of  the  bellmouth  at  the  inlet.  The  tunnel 
walls  were  set  with  a  total  of  3.8  cm  divergence  in  width  (51  cm  at  the 
entrance,  55  cm  at  the  exit).  Longitudinal  static  pressure  distributions 
were  measured  for  a  range  of  tunnel  speeds.  These  showed  no  systematic 
variations. 
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The  total  and  static  pressure  distributions  in  the  region  of  the  model 
locations  were  measured  in  the  second  phase  of  the  tunnel  calibration.  A 
rake  consisting  of  10  total  pressure  tubes  mounted  5  cm  apart  and  6  static 
pressure  tubes  mounted  5  cm  above  the  total  pressure  tubes  was  used.  The 
rake  data  were  taken  from  30  cm  above  the  model  location  to  40  cm  below 
the  model  location  at  5  cm  intervals.  The  resulting  pressure  data  are 
shown  in  Figure  19. 

The  tunnel  sidewall  static  pressure  was  measured  at  taps  on  each  sidewall 
approximately  10  cm  downstream  of  the  bellmouth  and  2.3  cm  above  the  tunnel 
floor  (on  the  model  chord  line).  This  sidewall  static  pressure  was  mea¬ 
sured  at  each  tunnel  set  condition  at  which  the  test  section  pressure  distri¬ 
butions  were  taken.  These  data  were  used  to  generate  the  tunnel  speed 
calibration  factor  which  relates  the  tunnel  dynamic  pressure  to  the  side- 
wall  static  pressure  drop  below  atmospheric  pressure.  This  calibration 

gaVfi  q  =  1.03  <Pq-Ps)  4-1 


POSITION  ACROSS  TUNNEL 


Figure  19.  Dynamic  Pressure  Distribution  At  Model 
Location 
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AIRFOIL  MODEL 


Model  Description 


The  airfoil  model  had  a  thickness  to  chord  ratio  of  0.20  and  was 
symmetrical  except  in  the  area  of  the  blowing  slot  and  upper  flap  contour. 
The  model  chord  was  20  cm  and  the  installed  model  span  was  51  cm.  A 
sketch  of  the  model  cross  section  is  shown  in  Figure  20  and  a  photograph 
is  shown  in  Figure  21. 


Figure  20.  Location  of  Pressure  Taps  on  Airfoil  Model 


Forward  of  the  maximum  thickness  ordinate  (X/C  =  0.4),  the  model 
contour  is  described  by  an  ellipse  which  gave  an  effective  leading  edge 
radius  of  approximately  5.57.  of  the  chord.  Aft  of  the  maximum  thickness, 
a  different  ellipse  was  used  to  describe  the  contour.  The  jet  blowing 
slot  was  located  on  the  upper  surface  at  X/C  =  0.9.  The  upper  nozzle  lip 
had  a  thickness  equal  0.0057.  of  the  model  chord.  The  107.  chord  flap  had 
a  circular  Coanda  turning  surface  with  a  radius  equal  to  1.4  cm  which  was 
tangent  to  the  flat  upper  surface  trailing  edge  at  X/C  =  0.93.  The  airfoil 
trailing  edge  was  slightly  blunt,  with  a  thickness  equal  0.017.  of  the  air¬ 
foil  chord.  The  lower  surface  trailing  edge  was  modified  slightly  from  the 
ellipse  to  give  a  sharp  edge  at  the  trailing  edge.  The  model  was  con¬ 
structed  of  aluminum.  Model  contours  in  the  center  5  cm  of  span  were  held 
to  +.005  cm. 

The  model  flap  was  attached  to  the  forward  portion  of  the  model  with 
one  of  five  separate  flap  brackets.  These  brackets  were  attached  inter¬ 
nally,  and  did  not  break  the  airfoil  mold  line.  The  flap  could  be 
mounted  at  angles  of  0°,  10°,  30°,  50°,  and  70°  relative  to  the  airfoil  chord 
plane.  The  flap  upper  surface  made  an  angle  of  20°  relative  to  the  chord. 
Thus,  the  flap  surface  and  hence  the  nominal  jet  deflection  angle  could 
be  set  to  20°,  30°,  50°,  70°,  and  90°  relative  to  the  airfoil  chord.  The 
effective  flap  pivot  point  is  at  X/C  =  .9  Y/C  =  0.05. 

The  total  span  of  the  constant  model  cross  section  was  61  cm.  The 
model  therefore  pierced  the  tunnel  walls  and  the  nozzle  gap  was  plugged 
to  the  correct  tunnel  span  of  51  cm.  The  nozzle  plugs  installed  at  each 
end  were  given  a  bellmouth  shape  in  the  planform  to  maintain  a  constant 
nozzle  total  pressure  distribution  up  to  the  wind  tunnel  walls. 

A  large  portion  of  the  model  interior  was  open  and  acted  as  the 
nozzle  feed  plenum.  The  flow  passage  from  the  plenum  to  the  nozzle 
convergence  was  designed  with  two  90°  turns  with  flow  areas  approximately 
twice  the  nozzle  area.  This  was  done  to  improve  the  spanwise  jet  total 
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pressure  distribution.  The  convergence  to  the  nozzle  exit  was  formed 
by  a  flat  upper  surface  and  the  circular  extension  of  the  Coanda  surface. 

The  flap  brackets  were  designed  such  that  the  nozzle  convergence  was  the 
same  for  all  flap  settings.  The  nozzle  exit  was  located  at  907.  chord 
and  designed  to  have  the  jet  exit  parallel  to  the  airfoil  chord. 

Ambient  temperature  compressed  air  was  supplied  to  the  model  through 
steel  endplates  to  which  were  welded  10  cm  pipe  fittings.  A  portion  of 
the  right  steel  endplate  could  be  removed  to  ease  flap  changes.  Sealing 
of  all  mating  surfaces  was  provided  by  0-rings. 

The  model  was  instrumented  with  a  total  of  82  pressure  taps. 

Seventy-four  taps  were  arranged  in  a  streamwise  row  on  the  airfoil  and 
flap  at  mid-span.  Of  these,  two  taps  were  located  on  the  blunt,  aft 
facing  surfaces,  the  nozzle  lip,  and  the  flap  trailing  edge.  An  addi¬ 
tional  eight  pressure  taps  were  arranged  in  two  spanwise  rows  to  monitor 
test  two-dimensionality.  Figure  20  shows  the  location  of  all  of  the 
pressure  taps. 

Model  Calibration 

The  model's  jet  nozzle  pressure  distribution  and  actual  jet  reaction 
force  and  angle  were  calibrated  prior  to  installation  in  the  two-dimensional 
tunnel. 

The  nozzle  was  plugged  to  the  correct  span  as  discussed  in  the  model 
description  section.  The  model  was  then  mounted  on  the  six  component 
pyramidal  balance  located  in  the  low  speed  wind  tunnel  at  Rockwell 
International's  Columbus  Plant.  The  air  supply  plumbing  to  the  model 
closely  matched  that  to  be  used  in  the  2-D  tunnel  installation.  A  venturi 
and  associated  instrumentation  were  installed  in  the  air  supply  line.  A 
pressure  tap  in  one  of  the  steel  model  end  fittings  was  used  as  the  cali¬ 
bration  reference  pressure.  One  objective  of  the  calibration  was  to 
generate  a  curve  relating  this  reference  pressure  to  the  nozzle  exit 
total  pressure. 
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Nozzle  exit  total  pressure  distributions  were  measured  at  each  flap 
angle  for  a  range  of  pressure  ratios  from  1.1  to  2.2.  These  distributions 
were  quite  uniform.  At  each  test  condition,  the  model  plenum  reference 
pressure  was  recorded.  These  data  were  then  used  throughout  the  test  to 
compute  exit  total  pressure. 

For  each  flap  angle,  a  set  of  force  data  were  taken  over  a  pressure 
ratio  range  of  1.1  to  2.1.  The  resultant  thrust  angle  was  computed  from 
the  model  lift  and  drag  forces  and  this  is  shown  in  Figure  22.  Model 
isentropic  thrust  was  calculated  from  the  venturi  measured  mass  flow  and 
the  nozzle  isentropic  velocity.  The  isentropic  velocity  was  computed 
assuming  an  expansion  from  the  nozzle  total  pressure  (obtained  from  the 
calibration  curve)  to  the  ambient  pressure.  Jet  temperature  was  measured 
at  the  venturi  as  experience  has  shown  this  to  be  within  a  few  degrees 


~  Flap  Upper  Surface  Angle  ~  Deg 
From  Horizontal 

Figure  22,  Comparison  of  Balance  and  Survey  Measured  Thrust  Angles 
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Tisentropic 


INSTRUMENTATION  AND  EQUIPMENT 


Survey  Probe 

Downstream  jet/wake  surveys  were  made  using  a  United  Sensor 
three-dimensional  directional  flow  probe,  type  DC- 125,  This  probe  has 
a  stepped  conical  sensing  head  with  one  pressure  tap  at  its  center  and 
four  taps  located  on  the  conical  surface.  This  probe  was  calibrated  as 
described  below  to  accurately  measure  total  and  static  pressure  and 
flow  angle  over  a  flow  angle  range  of  ±50  degrees  in  pitch  and  ±20  degrees 
in  yaw.  A  photograph  of  this  probe  is  shown  in  Figure  24. 


Figure  24.  Five  Port  Directional  Flow  Probe 


The  probe  was  calibrated  by  testing  it  in  a  jet  of  air  of  known 
total  and  static  pressure  and  flow  direction.  Data  were  taken  at  a 
matrix  of  probe  angles  of  +50  degrees  in  pitch  and  +20  degrees  in  yaw, 
in  5°  increments.  At  each  position,  the  five  separate  probe  pressures 
were  recorded.  These  data  were  repeated  in  calibration  jets  of  2.6, 
5.2,  and  7.75  cm  Hg  gage  total  pressure.  At  each  point,  the  following 
calibration  constants  were  calculated: 


KP  =  -  1  (+  Pitch,  i.e.,  PI  <  P3) 


KY  = 


P5-P3 

P3-P5 

P5-P1 

P5-P4 

P5-P2 

P2-P5 

P5-P4 


+  1  (-  Pitch,  i.e.,  PI  >  P3) 

-  1  (+  Yaw,  i.e.,  P4  >  P2) 

+  1  (-  Yaw,  i.e.,  P4  <  P2) 


vpT  =  P5-PT  (PT  =  jet  total  pressure 

PT-PS  PS  =  jet  static  pressure) 


KPS  = 


P5 


-P1+P2+P3+P4 


PT-PS 


These  data  are  then  used  to  generate  calibration  matrices  for  pitch 
angle  (function  of  KP  and  KY),  yaw  angle  (function  of  KP  and  KY),  total 
pressure  factor  (KPT  as  function  of  pitch  and  yaw  angle),  and  static 
pressure  factor  (KPS  as  function  of  pitch  and  yaw  angle). 


In  an  unknown  flow  field,  the  five  probe  pressures  ere  measured,  and 
used  to  compute  KP  and  KY,  These  values  are  used  as  the  independent 
variables  in  two  dimensional  interpolations  in  the  calibration  matrices 
to  give  the  pitch  and  yaw  angles  of  the  flow.  These  two  angles  are  then 
used  as  the  independent  variables  in  two  dimensional  interpolations  to 
give  the  flow's  total  and  static  pressure  factors,  KPT  and  KPS.  These 
factors  and  the  five  probe  pressures  are  used  to  compute  the  unknown 
flow's  total  and  static  pressure. 
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The  four  required  calibration  matrices  were  generated  using  the  data 
taken  in  the  7.75  cm  Hg  total  pressure  jet.  The  calibration  was  then 
checked  by  using  the  calibration  to  reduce  the  data  taken  at  the  other 
two  jet  pressures.  Flow  angles  were  generally  good  to  ±2  degrees. 

Total  static  pressures  were  within  +.05  cm  Hg. 

During  the  test,  the  3-D  probe  was  mounted  on  a  traversing  rig 
attached  to  the  tunnel  sidewall.  This  rig  consisted  of  two  support 
rails  and  a  motorized  screw  jack  which  could  position  the  probe  support 
arm  anywhere  along  76  cm  travel.  The  motor  was  under  either  manual  or 
computer  control  as  discussed  in  the  test  procedure  section.  The  tra¬ 
versing  rig  could  be  bolted  to  the  tunnel  wall  at  any  location.  The 
probe  could  be  positioned  at  any  spanwise  location  and  at  eight  discrete 
downstream  locations  for  any  rig  position. 

Instrumentation 

The  following  data  channels  were  used  during  the  2-D  wind  tunnel 
test! 

-  Wind  tunnel  sidewall  static  pressure 

-  Model  plenum  calibration  pressure 

-  Model  air  supply  venturi  upstream  pressure 

-  Model  air  supply  venturi  differential  pressure 

-  Model  air  supply  venturi  temperature 

-  Model  surface  pressures  -  5  module  scanivalve 

-  Jet/wake  3-D  survey  probe  -  5  low  range  and  5  high  range  transducers 

-  Jet/wake  3-D  survey  probe  traversing  rig  position 

-  Two-dimensional  wind  tunnel  ejector  supply  pressure 

-  Test  section  sidewall  boundary  layer  control  blower  plenum  pressure 

All  transducers  were  calibrated  prior  to  the  first  test  series  and 
were  recalibrated  prior  to  the  second  test  series. 
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All  instrumentation  output  signals  v  re  amplified  and  filtered  before 
going  to  the  multiplexed  analog  to  digital  converter  in  the  IBM  1800  data 
acquisition  computer.  Readings  were  then  converted  to  engineering  units 
and  preliminary  integrations  of  lift  (from  surface  pressures)  and  drag 
(from  both  wake  data  and  surface  pressures)  were  performed. 

The  model  surface  pressures  were  assigned  to  the  various  Scanivalve 
modules  on  the  basis  of  expected  maximum  pressures.  Each  Scanivalve  module 
had  the  transducer  reference  pressure  (the  tunnel  sidewall  pressure)  con¬ 
nected  to  the  first  and  last  ports.  Each  model  pressure  connected  to  the 
two  low  range  (±1  cm  Hg  max.  input)  modules  was  connected  to  two  consecu¬ 
tive  scanivalve  ports.  This  was  done  to  give  an  effective  settling  time 
of  twice  the  normal  settling  time  which  was  two  seconds. 

The  data  acquisition  program  was  written  to  read  every  other  port  on 
these  two  modules.  The  maximum  number  of  ports  read  was  set  by  the 
higher  range  modules  which  then  allowed  this  technique  to  be  used. 

The  3-D  wake  probe  was  required  to  measure  the  flow  properties  in  both 
the  low  total  pressure  tunnel  airflow  and  the  higher  total  pressure  jet 
sheet.  To  improve  the  accuracy  of  the  readings  in  the  low  pressure  flow, 
the  following  scheme  was  used.  Each  of  the  five  probe  pressures  were 
routed  to  two  transducers,  one  set  for  approximately  .775  cm  Hg  maximum 
input  (data  system  saturation  limit)  and  the  other  set  for  approximately 
26  cm  Hg  max  input.  A  pressure  relief  valve  was  installed  in  the  tube 
leading  to  the  low  pressure  transducer.  This  valve  was  set  to  shut  off 
this  pressure  line  at  about  1.5  cm  Hg,  well  below  the  2.5  cm  Hg  trans¬ 
ducer  physical  limit.  With  this  setup,  in  the  low  pressure  tunnel  flow, 
each  pressure  was  read  on  two  transducers,  and  the  data  reduction  program 
used  the  output  from  the  lower  range  transducer.  As  a  probe  entered  the 
jet,  the  low  range  transducer  data  channel  reached  saturation,  but  the 
relief  valve  protected  the  transducer  from  over  pressure  damage.  The 
data  reduction  program  recognized  this  saturation  value  and  instead  used 
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the  high  range  transducer  output  for  calculations.  As  the  probe  emerged 
from  the  jet,  the  relief  valve  reopened,  the  low  range  transducer  fell 
below  saturation,  and  its  more  accurate  output  was  again  used  for  that 
pressure.  The  over-lapping  ranges  of  transducer  electrical  saturation, 
relief  valve  operation,  and  transducer  physical  pressure  limits  allowed 
the  scheme  to  work  well  in  the  region  of  large  pressure  fluctuations  at 
the  edge  of  the  jet. 


TEST  PROCEDURES 


Data  Acquisition 

The  major  problem  with  two-dimensional  testing  is  the  effect  of  the 
sidewall  boundary  layer.  The  sectional  wing  lift  is  reduced  in  the  wall 
boundary  layer  which  gives  rise  to  three-dimensionality  and  the  associated 
induced  drag.  A  more  serious  effect  may  occur  when  the  wall  boundary  layer 
is  forced  to  separate  by  the  adverse  pressure  gradients  on  the  model.  The 
strength  of  these  adverse  gradients  is  significantly  increased  by  the 
large  lifts  generated  by  a  jet  flap  device. 

There  are  several  ways  to  mount  a  two-dimensional  model  on  a  balance, 
but  each  has  significant  problems  for  this  particular  test.  If  the  entire 
model  is  metric  (on  the  balance),  with  or  without  a  portion  of  the  sidewalls, 
the  measured  loads  will  include  the  non-uniform  load  near  the  walls  due  to 
the  boundary  layer.  In  this  case,  a  decision  must  also  be  made  whether  or 
not  to  attach  a  portion  of  the  walls  to  the  model.  In  the  first  case,  the 
tip  sealing  problem  is  alleviated,  but  the  forces  on  the  endplates  must  be 
accounted  for.  In  the  second  case,  a  force  free  seal  must  be  made  at  the 
model  tip/wall  juncture. 

A  second  approach  is  to  mount  only  a  center  span  portion  of  the  model 
on  the  balance.  This  eliminates  the  direct  force  effects  of  the  tip 
loading,  but  again  introduces  joints  and  seals,  near  the  region  of 
interest.  This  approach  is  substantially  complicated  on  a  jet  flap 
model  by  the  need  to  segment  the  nozzle. 

Any  balance  installation  on  a  jet  powered  model  faces  the  problem 
of  bridging  the  balance  with  the  air  supply  pipes.  This  problem  has 
been  successfully  solved  many  times  in  the  past.  But  each  installation 
is  generally  unique  and  there  is  usually  a  development  phase  before  a 
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useable  set-up  is  obtained.  In  any  case,  the  balance  must  be  calibrated 
with  the  air  supply  system  installed  and  pressure  effects  must  be  con¬ 
sidered.  All  these  effects  reduce  the  two-dimensionality  of  the  pressure 
distributions  and  measured  forces. 

In  view  of  the  above  complications,  the  decision  was  made  to  obtain 
lift  from  airfoil  surface  pressures  and  drag  from  wake  surveys. 

Integration  of  the  drag  component  of  the  surface  pressures  was  also  per¬ 
formed  as  a  secondary  drag  calculation.  Use  of  this  technique  also 
provides  information  on  the  qualitative  nature  of  the  flow. 

Each  run  during  these  tests  acquired  data  for  one  test  condition 
(angle  of  attack,  jet  deflection,  and  momentum  coefficient).  Model  angle 
of  attack  was  set  by  rotating  the  circular  wall  inserts.  The  wake  probe 
traversing  rig  was  positioned  such  that  the  probe  would  pass  the  desired 
distance  downstream  of  the  model  (approximately  15  cm).  The  rig  tilt 
angle  was  set  such  that  the  expected  flow  angles  would  be  within  the 
probe's  acceptance  limits  (±50  degrees  from  the  probe  axis).  The  desired 
tunnel  speed  was  set  by  adjusting  a  pressure  regulating  valve  while  moni¬ 
toring  the  sidewall  static  pressure  displayed  on  a  digital  voltmeter  (DVM) 
Model  thrust  was  set  by  operating  a  pressure  control  valve  to  achieve  the 
precomputed  model  plenum  reference  pressure  also  displayed  on  a  DVM. 
Sidewall  BLC  operating  pressure  was  set  in  the  same  way  as  the  model 
pressure.  Proper  set  pressures  were  determined  by  observing  tufts 
attached  on  the  sidewall  and  model  upper  surface  and  by  monitoring  the 
model  spanwise  row  of  pressure  taps. 

With  the  desired  test  conditions  set,  data  acquisition  was  under 
computer  control.  Model  and  tunnel  operating  conditions  were  read 
followed  by  the  scanivalve  data.  The  computer  then  commanded  the  tra¬ 
versing  rig  travel  based  on  a  predetermined  number  of  datapoints  and 
desired  distance  between  points.  Different  numbers  of  points  and  point 
spacing  could  be  requested  in  each  of  three  survey  regions  (nominally 
above  the  jet,  in  the  jet,  and  below  the  jet).  After  each  step  of  probe 
travel,  a  pressure  setting  time  of  two  seconds  was  allowed  before  data 
were  acquired. 
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For  each  channel  of  data  acquired,  the  computer  took  ten  readings 
over  a  V2  second  period.  These  readings  were  averaged  and  this  average 
was  stored  for  later  conversion  to  engineering  units. 

Each  run  required  approximately  15  minutes  and  required  about  190  wake 
probe  positions.  Test  conditions  were  held  constant  during  this  time  by 
two  technicians.  The  test  engineer  was  able  to  monitor  these  conditions 
and  progress  of  the  test  from  the  data  acquisition  system. 

Two  test  series  were  run  several  months  apart.  The  objective  of  the 
first  series  was  to  obtain  lift  and  drag  data  for  jet  deflection  angles 
of  30°,  50° i  70°,  90°,  and  momentum  coefficients  between  0.5  and  24.0. 

The  objective  of  the  second  test  was  to  survey  the  jet  at  five  down¬ 
stream  locations  for  a  range  of  test  conditions  to  quantify  the  jet 
thrust  and  mass  flow  changes  as  the  jet  develops  downstream. 

Data  Reduction  Procedures 

The  following  section  discusses  the  data  reduction  procedures  for 
the  wind  tunnel  and  model  operating  conditions,  airfoil  surface  pressure 
lift  and  drag  integrations,  wake  drag  integration,  and  jet  profile  thrust 
and  mass  flow  integrations. 

The  wind  tunnel  operating  conditions  are  computed  from  the  measured 
sidewall  static  pressure  and  ambient  pressure  and  temperature.  Because 
of  the  suck  down  type  tunnel  design,  the  test  section  total  pressure  is 
equal  to  the  ambient  pressure;  this  was  confirmed  by  the  test  section 
calibration. 

^Test  Section  P ambient 

This  calibration  provided  the  following  relationship  between  the  sidewall 
static  pressure  and  the  calibrated  test  section  uncorrected  dynamic 
pressure. 

^test  section  uncorrected  =  x  ^Psidewall  ”  Pambient^ 
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A  solid  blockage  correction  was  included  to  give  the  corrected  test 

section  dynamic  pressure. 

q„  .  .  .  =  q  (1.0  +  .0006  )2 

test  section  ^uncorrected 

Test  section  static  pressure  is  then: 

Ptest  section  “  ’  ^TS 


Test  section  Mach  number 

“ts  * 


'P-A  .28571 

J/ 


x  5 


Model  Reynold's  number  (based  on  model  chord): 


Rn  =  106  (1 


.798)  ,TS  [TTS  *  ^  j 


4-2 


4-3 


Where  TTg  is  ambient  temperature. 


The  model  jet  isentropic  thrust  was  computed  from  venturi  measured 
mass  flow  and  jet  isentropic  velocity.  The  venturi  data  (PUpstream» 
differential  *  and  Temperature)  were  reduced  to  a  mass  flow  using 
standard  computations.  Isentropic  velocity  was  computed  as  follows: 


The  nozzle  exit  total  pressure  was  obtained  from  the  model  plenum 
calibration  pressure  and  the  calibration.  The  jet  temperature  was  assumed 
to  be  equal  to  the  measured  venturi  temperature  as  discussed  earlier.  The 
isentropic  velocity  is  then  calculated  assuming  isentropic  expansion  from 
the  nozzle  total  pressure  to  the  test  section  static  pressure. 


isen 


dest  Section 
^nozzle 


.28571 


4-4 


Where  7  is  the  ratio  of  specific  heats  and  R  is  the  universal  gas 


constant. 


Each  of  the  wake/jet  sheet  surveys  from  both  test  series  has  been 
integrated  to  give  an  airfoil  drag  coefficient.  This  integration  was 
performed  using  the  method  of  Betz,  (Schlichting,  ref  28),  modified  to 
include  the  effect  of  the  large  flow  deflection  angle.  Direct  application 
of  the  Betz  method  of  wake  survey  integration  does  not  give  the  correct 
value  of  drag,  if  the  survey  is  made  before  the  jet  has  completely  turned 
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to  the  free  stream  direction.  In  that  case,  there  is  a  static  pressure 
difference  which  makes  a  contribution  to  the  horizontal  force.  This  con¬ 
tribution  must  be  equal  to  the  difference  between  the  horizontal  component 
of  the  local  jet  thrust  and  the  thrust  of  the  completely  deflected  jet. 
Since  the  jet  approaches  the  free  stream  direction  asymptotically,  the 
actual  thrust  can,  in  principle,  only  be  determined  from  a  jet  survey  at 
downstream  infinity,  or  from  a  survey  which  extends  to  infinity  above  and 
below  the  jet.  However,  if  the  jet  is  nearly  parallel  to  the  free  stream, 
the  tangential  component  of  the  jet  thrust  is  equal  to  the  value  recovered 
at  infinity,  within  experimental  error.  This  will  be  seen  in  the  analytic 
results.  Thus,  the  measured  tangential  thrust  was  used  to  determine  the 
thrust  recovery. 
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V.  RESULTS  AND  DISCUSSION 

SURFACE  PRESSURE  DISTRIBUTIONS 

The  character  of  the  measured  surface  pressure  distributions  changes 
as  the  jet  thrust  coefficient  and  deflection  angle  are  increased.  For 
low  thrust  coefficients  (Cp  <  1)  and  small  deflection  angles  (9  =  30°), 
the  pressure  distributions  have  a  saddle  shape,  typical  of  that  predicted 
by  classical,  inviscid  jet  flap  theory.  A  representative  case  is  shown  in 
Figure  25.  The  saddle  shape  is  due  to  the  suctions  peaks  at  the  leading 


X/C 

Figure  25.  Surface  Pressure  Distribution  for  Attached  Flow,  Cy  *  2,  0O  * 
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edge  of  the  airfoil  and  the  knee  of  the  flap.  The  pressure  is  further 
reduced  on  the  upper  surface  of  the  flap  by  the  suction  which  deflects 
the  jet  over  the  Coanda  surface.  Since  the  radial  pressure  gradient 
balances  the  inertia  of  the  jet,  the  reduction  in  flap  surface  pressure 
is  given  by 


ACp  =  Cpc/R  5-1 

The  sudden  pressure  rise  at  the  point  where  the  surface  of  the  flap  changes 
from  a  circular  arc  to  a  straight  section  causes  a  separation  bubble  to 
develop;  the  suction  peak  near  the  nozzle  exit  is  due  to  this  bubble.  The 
smaller  rise  at  the  trailing  edge  of  the  flap  is  due  to  the  jet  flap  effect: 
The  jet  leaves  the  airfoil  along  a  tangent  to  the  upper  surface  of  the  flap, 
and  therefore  deflects  the  trailing  edge  streamline  relative  to  the  mean 
line  of  the  airfoil.  This  produces  a  suction  peak  similar  to  the  one  which 
occurs  at  the  knee  of  a  flap. 

For  moderate  values  of  the  thrust  coefficient  (2  <  Cy  <  8)  and  deflection 
angles  (6  *  30°,  50°),  a  separation  bubble  appears  on  the  upper  surface  of 
the  airfoil.  The  flow  separates  from  the  leading  edge  of  the  airfoil  and 
reattaches  at  some  point  on  the  upper  surface.  In  this  case,  the  pressure 
is  nearly  constant  over  the  front  of  the  airfoil,  but  returns  to  the  ideal 
saddle  shape  over  the  rear  section.  A  typical  example  is  shown  in  Figure 
26. 


For  large  values  of  the  thrust  coefficient  (Cy  >  8)  and  large  deflection 
angles  (9  >  50°)  the  separation  bubble  does  not  reattach  to  the  surface  of 
the  airfoil.  In  this  case,  the  pressure  distribution  is  flat  or  even 
slightly  arched,  as  seen  in  Figure  27.  Since  the  combined  effect  of  the 
jet  and  flap  accelerates  the  flow  over  the  unblown  section  at  the  rear  of 
the  airfoil,  trailing  edge  separation  did  not  occur  for  any  case  tested. 

The  complete  set  of  surface  pressure  distributions  may  be  found  in 
Appendix  I. 
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WAKE  PRESSURE  DISTRIBUTIONS 

Distributions  of  the  total  and  static  pressure,  as  well  as  the  flow 
angularity,  were  obtained  in  the  near  wake.  The  trend  of  the  wake  data 
was  consistent  with  the  surface  pressure  distributions.  A  representative 
set  of  pressure  profiles  for  the  case  in  which  the  flow  did  not  separate 
from  the  airfoil  is  shown  in  Figure  28.  The  total  pressure  profile  has 
the  classical  bell  shape,  and  the  static  pressure  is  less  above  the  jet 
than  below  it.  In  Spence's  (ref  10)  theory,  this  pressure  jump  is 
represented  by  the  vortex  sheet.  The  reduction  in  static  pressure  within 
the  jet  is  due  to  the  jet  turbulence,  which  causes  an  actual  reduction 
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Figure  28.  Total  and  Static  Pressure  Surveys  in  the  Jet  Wake  for  Cu  =>  1 
and  0  =  50° 
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case  with  Cy  =  1.  As  might  be  expected,  the  penetration  of  the  jet 
increases  with  the  initial  deflection  angle.  The  effect  of  the  jet 
momentum  coefficient  on  the  trajectory  is  shown  in  Figure  30  for  two 
values  of  the  initial  deflection  angle.  As  the  blowing  coefficient  is 
increased,  the  jet  tends  to  straighten  out.  In  the  limit  as  C)j  -*■  °°,  thi 
jet  continues  along  the  initial  trajectory  and  the  horizontal  force  on 
the  airfoil  becomes  TcosQ.  These  wake  survey  data  are  contained  in 
Appendix  III. 


X/C 

Figure  30.  Variation  of  Jet  Shape  with  Jet  Momentum,  Cu 


LIFT  AND  THRUST  INTEGRALS 


The  measured  surface  pressure  distributions  were  integrated  to  obtain 
the  lift  on  the  airfoil.  In  Figure  31  the  total  lift  coefficient  is  shown 
as  a  function  of  the  jet  thrust  coefficient.  The  total  lift  on  the  airfoil 
is  equal  to  the  sum  of  the  circulation  lift  and  the  vertical  reaction  to 
the  jet  thrust,  which  appears  as  a  suction  on  the  flap  surface.  For  the 
smaller  jet  deflection  angles  (30°,  50°)  the  data  are  well  represented  by 
the  theoretical  results  obtained  by  Spence  (ref  10).  However,  for  the 
larger  deflection  angles,  there  is  a  loss  of  circulation  lift  due  to 
separation  from  the  airfoil.  Impingement  of  the  jet  on  the  tunnel  floor 
causes  the  lift  to  increase  again  at  large  values  of  and  9,  but  this  is 
a  blockage  effect  and  not  representative  of  the  force  that  would  be  induced 
by  an  unconfined  jet. 


Figure  31.  Lift  Coefficient  as  a  Function  of  Momentum  Coefficient 


In  order  to  provide  a  check  on  the  experimental  methods,  the  thrust 
recovery  was  calculated  from  both  the  wake  survey  data  and  the  airfoil 
surface  pressures.  As  previously  explained,  the  method  of  Betz  (Schlichting , 
ref  28)  was  used  to  evaluate  the  excess  momentum  integral  from  the  wake 
pressure  distributions.  The  net  thrust  on  the  airfoil  is  then  obtained 
by  adding  the  source  thrust  of  the  initial  jet  mass  to  the  momentum 
integral.  Since  the  thrust  recovery  is  defined  relative  to  the  viscous 
drag  of  the  unblown  airfoil,  this  baseline  drag  is  added  to  the  net  thrust 
of  the  jet  flap.  Thus,  the  loss  of  recovery  is  taken  to  be  the  increase 
in  the  profile  drag  of  the  unblown  airfoil  due  to  jet  entrainment  and 
flow  separation  induced  by  the  jet.  The  wake  recovery  factor  is  there¬ 
fore  given  by 


Tx  +  Tq+|  Dv 

rw  *  - f -  5-2 

in  which  Tx  ic  the  excess  momentum  integral,  Tq  is  the  source  thrust,  and 
Dv  is  the  viscous  drag. 

The  viscous  drag  of  the  airfoil  was  estimated  by  computing  the  pressure 
and  velocity  distributions  on  the  unblown  airfoil.  The  actual  airfoil 
coordinates  and  test  Reynolds  number  were  used.  Variations  in  lift  were 
obtained  by  changing  the  angle  of  attack.  In  Figure  32  the  predictions 
of  this  program  are  compared  to  measurements  for  10°  angle  of  attack.  The 
excellent  agreement  of  the  measured  and  predicted  surface  pressure  distribu¬ 
tions  is  an  indication  that  good  two  dimensionality  was  achieved  in  the 
wi:  \  tunnel  test.  The  agreement  between  the  lift  coefficients  is  also 
good.  The  measured  axial  force  coefficient  is  within  10%  of  the  computed 
value;  since  integration  of  the  surface  pressure  data  to  obtain  the  axial 
force  is  difficult,  this  is  considered  good  agreement.  The  computed 
viscous  drag  coefficient  of  the  airfoil  is  described  by  the  function 

Cv  =  .007Cl2  +  .022  5-3 
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Figure  32.  Comparison  of  Measured  and  Predicted  Surface  Pressure 
Distributions 

in  which  Cl  is  the  circulation  lift  coefficient.  The  viscous  drag 
includes  both  the  skin  friction  drag  and  the  pressure  drag  due  to  the 
displacement  effect  of  the  boundary  layer. 

The  drag  obtained  by  integrating  the  airfoil  surface  pressures  also 
includes  the  pressure  drag  due  to  viscosity,  but  not  the  skin  friction 
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drag.  Thus,  the  pressure  drag  of  the  unblown  airfoil  is  added  to  the 
measured  surface  pressure  drag.  In  addition,  however,  an  induced  drag 
correction  must  be  made  because  end  effects  introduce  some  three 
dimensionality  in  the  surface  pressure  distributions.  The  airfoil 
recovery  factor  is  therefore  given  by 


ra  -  1  - 


i  Ds  +  Dp  +  Pi 

;  T 


5-4 


in  which  Ds  is  the  integral  of  the  surface  pressures.  Dp  is  the  computed 
viscous  pressure  drag  of  the  unblown  airfoil,  and  Di  is  the  induced  drag 
correction. 


Foley  (ref  2)  estimated  the  induced  drag  by  extrapolating  the  linear 
section  of  the  drag  polar  of  the  unblown  airfoil  to  the  large  lift  coef¬ 
ficients  of  the  jet  flapped  airfoil.  This  is  equivalent  to  assuming  that 
the  lift  and  jet  thrust  are  elliptically  distributed  along  the  span.  A 
similar  approach  was  taken  in  this  study,  but  the  large  values  of  the 
thrust  coefficient  considered  here  required  using  the  expression  for  the 
induced  drag  of  a  jet  flapped  wing  devised  by  Maskel  and  Spence  (ref  7), 

CL2 

°Di  =  nARe  +  2Cp  5-5 

In  effect,  Foley  (ref  2)  neglected  the  thrust  coefficient  term,  but  over 
the  range  of  values  he  examined  (Cy  <  1),  the  error  is  small.  In  this 
study,  the  effective  aspect  ratio  of  the  airfoil  was  determined  by  sub¬ 
tracting  the  computed  viscous  drag  from  the  measured  total  drag  of  the 
unblown  airfoil.  The  measured  drag  polar  is  shown  in  Figure  33.  The 
total  drag  is  described  by  the  function 

CD  =  .011CL2  +  .022  5-6 

The  effective  aspect  ratio  of  the  airfoil  is  then  obtained  from  the 
relation 
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Figure  33.  Drag  Polar  of  the  Unblown  Airfoil 

“*  •  .(C„C-  C)  5'7 

This  gives  ARe  =  77,  which  was  then  assumed  to  be  independent  of  variations 
in  Cy  and  90. 

The  values  of  the  thrust  recovery  factor  obtained  by  both  methods  are 
compared  in  Figure  34.  There  is  surprisingly  good  agreement  between  the 
methods,  although  the  surface  pressure  data  gives  consistently  higher 
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Figure  34.  Variation  of  the  Thrust  Recovery  Factor 

recovery.  Based  on  our  prior  experience  with  wake  and  surface  pressure 
integrations,  it  is  our  feeling  that  the  wake  survey  data  is  more  reliable 
However,  the  differences  observed  are  within  the  16%  range  of  experimental 
error  associated  with  previous  measurements  of  thrust  recovery  (Leaman 
and  Plotkin,  ref  29).  The  increase  in  the  recovery  factor  observed  for 
large  values  of  Cy  and  90  is  due  to  jet  impingement  on  the  tunnel  floor, 
and  is  not  representative  of  the  thrust  recovery  of  an  unconfined  jet. 

ANALYSIS  OF  RESULTS 

In  general,  the  thrust  recovery  decreases  as  the  jet  thrust  coefficient 
and  initial  deflection  angle  are  increased.  However,  the  variation  is  not 


continuous.  For  every  value  of  the  deflection  angle  there  is  a  discon¬ 
tinuous  change  in  the  recovery  for  a  value  of  the  thrust  coefficient  of 
approximately  Cy  =  2.  During  the  test,  a  distinct  increase  in  noise 
level  could  be  heard  when  the  discontinuity  occurred.  Comparison  of  the 
measured  airfoil  surface  pressure  distributions  just  before  and  after  the 
discontinuity  reveals  that  it  is  due  to  the  bursting  of  the  leading  edge 
separation  bubble.  This  is  shown  in  Figure  35  for  the  90  =»  50°  case. 
There  is  a  change  between  a  short  separation  bubble,  which  gives  larger 
recovery  and  a  long  bubble  which  gives  less  recovery.  Although  the 
thrust  recovery  continues  to  decrease  as  the  long  bubble  grows,  there  is 
no  change  when  the  bubble  ultimately  fails  to  reattach  to  the  surface  of 


X 


Figure  35.  Comparison  of  Measured  and  Ideal  Airfoil  Pressure  Distributions 


Che  airfoil  at  all.  Because  earlier  studies  were  conducted  at  small 
values  of  Cy,  this  discontinuity  has  not  been  previously  known. 


The  bursting  of  the  bubble  depends  on  the  Reynolds  number  and  Mach 
number  of  the  flow,  and  the  leading  edge  radius  of  the  airfoil,  so  that 
these  trends  cannot  be  generally  applied.  However,  the  value  of  the  thrust 
recovery  factor  is  well  predicted  by  the  jet  flap  analysis  up  to  the  point 
of  bursting.  Thus,  it  can  be  concluded  that  the  loss  of  recovery  up  to 
the  burst  point  is  due  to  jet  entrainment  drag,  and  that  the  short  bubble 
does  not  significantly  affect  either  the  lift  or  thrust  integrals. 


In  Figure  36  the  computed  surface  pressure  distribution  on  the  unblown 


undeflected  jet,  for  Cp  =  1.  The  effect  of  the  jet  is  to  increase  the 
base  drag  on  the  airfoil.  Integration  of  the  surface  pressures  gives 
r  =  0.050  for  the  thrust  recovery.  The  corresponding  variation  in  the 
excess  momentum  of  the  jet  is  shown  in  Figure  37.  By  the  last  station 
shown,  the  jet  thrust  has  decreased  to  give  r  =  0.051. 

For  small  values  of  the  jet  thrust  coefficient,  the  thrust  recovery 
actually  increases  slightly  as  the  jet  deflection  angle  is  increased  from 
30°  to  50°.  Similar  behavior  was  seen  in  Foley's  (ref  2)  data,  and 
predicted  by  the  jet  flap  analysis.  Because  the  static  pressure  on  the 
upper  surface  of  the  jet  becomes  more  negative  as  the  jet  is  deflected, 
there  may  actually  be  some  thrust  augmentation  in  this  case.  However, 
the  effect  is  small  and  not  sufficient  to  overcome  the  entrainment  drag 
which  occurs  farther  downstream. 
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Figure  37.  Streamwise  Variation  of  Excess  Momentum  in  Jet 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  thrust  recovery  is  generated  by  equal  but  opposite  pressure  forces 
induced  on  the  airfoil  and  the  jet  flap.  If  the  mixing  of  the  jet  is 
neglected,  the  jet  is  a  streamline  of  the  flow,  and  the  net  pressure  force 
is  everywhere  normal  to  the  jet.  In  this  case,  the  thrust  of  the  jet  is 
conserved  and  ultimately  recovered  as  a  horizontal  force  in  the  far  wake. 

This  gives  complete  thrust  recovery.  Real  fluid  effects  reduce  the  actual 
recovery  both  by  inducing  a  jet  drag  on  the  airfoil  and  by  causing  separation 
of  flow  from  the  airfoil. 

The  thrust  recovery  is  a  function  of  both  the  thrust  coefficient  and  jet 
deflection  angle.  In  general,  the  recovery  is  nearly  complete  for  small 
values  of  the  thrust  coefficient  (Cy  <  1),  but  decreases  to  the  horizontal 
component  of  the  jet  reaction  force  at  large  values  of  thrust  (Cy  >  15). 

The  effect  of  jet  .mixing  is  to  reduce  the  thrust  recovery  about  10%  at 
small  values  of  the  thrust  coefficient  and  deflect*'  angle.  As  these 
parameters  are  increased,  flow  separation  causes  a  further  loss  in  recovery, 
t  For  a  fixed  deflection  angle,  there  is  a  discontinuous  change  in  the  recovery 

factor  at  Cy  -  2,  as  the  character  of  the  separation  changes  from  a  "short" 
bubble  which  reattaches  to  the  airfoil  near  the  leading  edge,  to  a  "long" 
bubble  which  reattaches  near  the  trailing  edge.  The  flow  is  bistable  in  the 
transition  region.  At  higher  values  of  Cy,  the  bubble  does  not  reattach  to 
the  airfoil  at  all,  but  is  entrained  into  the  wake. 

A  vlscid/inviscid  interaction  analysis  has  been  developed  to  compute 
the  contribution  of  the  jet  drag  to  the  reduction  in  thrust  recovery. 

It  was  found  that  the  entrainment  drag  accounts  for  all  the  variation  in 
recovery  before  bursting  occurs.  Thus,  such  an  analysis  can  be  used  to 
predict  the  thrust  recovery  of  various  jet  flap  concepts  if  the  flow  is 
prevented  from  separating  by  the  use  of  slats,  etc. 

Because  flow  separation  plays  such  an  important  part  in  determining  the 
thrust  recovery  at  large  values  of  the  thrust  coefficient,  further  develop¬ 
ment  of  jet  flap  theory  should  focus  on  predicting  the  effects  of  separation. 
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0,  Theta  *=  0°,  Alpha  -  -5°  (Test  272,  Run  3) 
7.28  psf 


0,  Theta  =  0°,  Alpha  -  15°  (Test  272,  Run  8) 
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Figure  1-6.  C  -  0.53,  Theta  -  30°,  Alpha  -  0°  (Test  272,  Run  9) 
Q  ■  6.91  psf 
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1.07,  Theta  =  30°,  Alpha  =  0°  (Test  272,  Run  10) 
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Figure  1-11.  Cy  =  6.40,  Theta  =  30°,  Alpha  =  0°  (Test  272,  Run  15) 
Q  =  4.04  psf 
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Figure  1-16.  ■  2.05,  Theta  ■  50°,  Alpha  *  0°  (Test  272,  Run  24) 

Q  -  7.02  psf 


I  w  -4)  -©  m?r  ^lao 

•  *  '  *  t  ,  i  . 


Figure  1-21.  C  -  8.17,  Theta  =  50°,  Alpha  *  0°  (Test  272,  Run  25) 
Q  -  3.02  psf 
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Figure  1-22.  C  -  16.22,  Theta  =  50°,  Alpha  -  0°  (Test  272,  Run  26; 
Q  ■  1.51  psf 


- 9<S&— — q  — o~0(P1q 


70°,  Alpha  =  0°  (Test  272,  Run  43) 
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Figure  1-36.  Cy  *  0.97,  Theta  =  90°,  Alpha 
Q  =■  6.95  psf 
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WAKE  SURVEY  DATA 

TEST  27Z  RUN  _A_ 
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INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 
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Figure  II-4. 
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TEST  2J1  RUN  _9_ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  J£_  C  MU 


»«•!*  M.K  14.  M  14. 3S  14.4* 


PRESSURE  ~  PS I A 


Xss/ chond  —  ■  3  73  3  Yas/chond  —  *  ^ 
BETA  =  2^03 


Figure  II-7. 


ALPHA  ~  DEGREES 


141 


WAKE  SURVEY  DATA 

TEST  RUN  / o 

ANGLE  OF  ATTACK  _Q__  NOMINAL  JET  ANGLE  2£_  C  MU  L°7 


Xse/chord  *  >3?*/  Yss/chord  =  »  S*i 3 3 
BETA  =  30-33 


142 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 
TEST  RUN  J_L 

ANGLE  OF  ATTACK  _2_  NOMINAL  JET  ANGLE  C  MU 


Xso/chord  =  >3581  Yea/ chard  ~  llQJJo 
BETA  =  2o.&  3 
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WAKE  SURVEY  DATA 
TEST  .272  RUN  /  3 


PRESSURE  ~  PSIA  j 

.  i 

Xoe/chord  58  •  35~^/  Yee/chord  53  -/Of 6  | 


BETA  =  2Q-B3 


WAKE  SURVEY  DATA 


TEST  ZJ&  RUN  >  y 

ANGLE  OF  ATTACK  O  NOMINAL  JET  ANGLE  3 O  r  MU 


PRESSURE 

Xso/chor'd  =  -it  8 / 

f  Yoo/chond  =»  -/032 

BETA 

-  ZO-B3 

•  |  »  »  i  ■  |  i  i  •  i  |  »  •  »  i  |  i 

•  5  J*  1* 


1*11111 


Figure  11-11. 


M  H 

ALPHA  ~  DEGREES 
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INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 
TEST  Hi  RUN  JA. 

ANGLE  OF  ATTACK  _Q_  NOMINAL  JET  ANGLE  JA  C  MU  8.V398 


- 

- 

- 

- 

\ 

- 

J _ 

1 

I 

1 

I 

14.*  14.1  14.4  14.*  14.*  IS.* 


PRESSURE  ~  PSIA 

Xse/ chord  23  <  35* 6  9  Yoo/chord  *  >  0  9 BS~ 

BETA  =  gg-83 
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WAKE  SURVEY  DATA 

TEST  ZTZ  RUN  J  7 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  _2 0_  C  MU  8.60/? 


14.2  14.2  14.4  14.2  14.2  1S.2  12.2 


PRESSURE  "  PS I A 

Xse/choi  d  53  »  3&A 9  Yee/chord  —  ./OOQ 

BETA  =  2083 


I*  -2  •  2  12  12  22  22  M 


ALPHA  ~  DEGREES 


Figure  11-14. 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  ZJ2,  RUN  JJL 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET  ANGLE  £Q_  C  MU  /.  Q/3 


Xee/chord  58  -Z2SZ  Ye®/ chord  —  .QS/  3 
BETA  =  30 
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WAKE  SURVEY  DATA 

TEST  HI  RUN  ZL 

ANGLE  OF  ATTACK  _jO_  NOMINAL  JET  ANGLE  _££.  C  MU  2.0 vr 


>4.2  14.3  14.4  14. S  14. t  14.7  14. B 


PRESSURE  ~  PS I A 

Xee/chord  -  •  ZOV/  Ysa/^hord  *  “-Q/ 66 
BETA  -  So 


■«  •  S  IB  IS  IB  IS 


Figure  11-18. 


ALPHA  ~  DEGREES 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  RUN  _2_V 

ANGLE  OF  ATTACK  _C_  NOMINAL  JET  ANGLE  _££_  C  MU  £ 


Xae/chord  =  Yes/ chord  =  •  °3&9 

BETA  »  3g. 


154 


.  os^rh 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  372  RUN  AL 

ANGLE  OF  ATTACK  __Q_  NOMINAL  JET  ANGLE  J[£_  C  MU  306VT 


PRESSURE  ~  PS I A 

Xso/ohord  -  . 2</P9  Yse/chord  =  -0386 
BETA  =  30 
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INCHES 


WAKE  SURVEY  DATA 

TEST  2J2  RUN  2Z_ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  C  MU 


WAKE  SURVEY  DATA 

TEST  2JZ  RUN 

ANGLE  OF  ATTACK  _2_  NOMINAL  JET  ANGLE  _££  C  MU  g. !Z°Q 


Xaa/chord  =  •  2V  Ob  Yaa/chord  =  .0379 
BETA  -  30 
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INCHES 


WAKE  SURVEY  DATA 

TEST  il£  RUN  2£_ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  £0_  C  MU  /£JL2'3 


X««/chord  “  >  3  '5Q£>  Yse/chord  -  ■  0&Q5 
BETA  =  VO 


Figure  11-24. 


INCHES 


WAKE  SURVEY  DATA 


TEST  272  RUN  _2J_ 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET  ANGLE  JjO_  C  MU  22-33  V 6 


Xaa/chord  -  -  B5Q9  Yaa/chord  “  -02/  2 


BETA  »  30 


x 

M 
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WAKE  SURVEY  DATA 
TEST  221  RUN  30 

ANGLE  OF  ATTACK  JL_  NOMINAL  JET  ANGLE  !£_  C  MU 


BETA  =  «f.  2- 


ALPHA  "  DEGREES 

Figure  11-26. 


3 


WAKE  . SURVEY  DATA 

TEST  ill  RUN  _V3. 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  70_  C  MU  '■ 0 v^6 


Xso/chord  *»  >192- 1>  Y»»/chord  -  -.03 S~3 
BETA  « 


161 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  273  RUN  3/ 

ANGLE  OF  ATTACK  _JL  NOMINAL  JET  ANGLE  LQ_  C  MU  *  0**3 

12. S  !  I  I  I  I - - ! - 1 - 1 - 1 - 1 


WAKE  SURVEY  DATA 


Xee/chord  =  •l*f7Z  Yea/chord  »  ~-0 / 8 
BETA  =  V5-.Z 


Figure  11-33. 


DEGREES 


167 


i 

f- 

k 


WAKE  SURVEY  DAJA. 

TEST  2J2  RUN  _1£ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  70_  C  MU  28. 9876 


Xso/chord  **  •  / 2  2  Ysa/chond  »  “>033/ 
SETA  -  *£2 


170 


WAKE  SURVEY  DATA 

TEST  272  RUN  Vv 


BETA  »  Z 


13.5 


1*.» 


CO 

u; 

-i. 

u 

z 


T.S 


i 

JSI  5* 


3.5 


ALPHA 


DEGREES 


Figure  11-37. 
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INCHES  2  ~  INCHES 


WAKE. SURVEY  DATA 

TEST  ZUZ  RUN  JUL 


ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  OJL  C  MU  -520/ 


Xa«/ehor~d.  —  >J7 5* 9  Y©a/chord  —  •  O  V/6 
BETA  *  Vi'g 


172 


WAKE  SURVEY  DATA 
TEST  ill  RUN  X&. 


Figure  11-41. 


ALPHA  ~  DEGREES 


175 


WAKE  SURVEY  DATA 

TEST  2JZ  RUN  Jt± 

ANGLE  OF  ATTACK  1 Q_  NOMINAL  JET  ANGLE  ±Q_  C  MU  //  *979 


X®e/chord  -  >99  ty  Yae/chord  =  •P.5‘73 


| 


Figure  11-43. 


Figure  11-44. 
178 


ALPHA  ~  DEGREES 


INCHES  2  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  %JJi  RUN 

ANGLE  OF  ATTACK"  O  NOMINAL  JET  ANGLE  JLSL  C  MU  Z£Jf* *  6 


PRESSURE 

Xee/  chord  —  « S3 8 V  Ys»/ chord 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  g&l  RUN  V 2 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET  ANGLE  _3 C  MU  '.0*09 


- 

L_ 

- 

- 

L 

-4 

1 

] 

14. 

\ 

9  14. 

1 

4  14. 

- 1 - 

6  14. 

I 

t  15. 

i 

I  15. 

J 

2  15.4 

PRESSURE  ~  PS I A 


Xss/chord  =  Yss/chord  = 

beta  =  / O 


182 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA. 

TEST  231  RUN  H  3 

ANGLE  OF  ATTACK  _Q_  NOMINAL  JET  ANGLE  _3j2_  C  MU  1 


PRESSURE  ~  PSIA 

Xee/chord  —  • V/ 3  Yss/ chord  -  •  7 38 

BETA  =  / Q 


Figure  III-2. 

183 


INCHES 


WAKE  SURVEY  DATA 

TEST  18/  RUN  %£_ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  JO_  C  MU  /os'9 


Xaa/chord  -  ■ 783  Yea/ chord  “  -^73 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  ISZ  RUN  V7 

ANGLE  OF  ATTACK  _Q_  NOMINAL  JET  ANGLE  C  MU 


1«.S  IS.*  SS.S  IS.*  IC.S 


PRESSURE  ~  PS I A 

Xee/chord  =  _« Z  9  Yea/chord  -  •  76 

BETA  =  / O 


Figure  Hl-6. 


ALPHA  ~  DEGREES 


INCHES 


WAKE  SURVEY  DATA 

TEST  232  RUN  _v£ 

ANGLE  OF  ATTACK  JL.  NOMINAL  JET  ANGLE  2Q_  C  MU  2-086,2 


- 

5 

j  _ 

■ 

^ — r-i — | — t — t — - — i — 1 — i — ■ — t — t — I 

14.4  14. S  14.1  IS.*  15.2  IS.  4  1S.6 


PRESSURE  ~  PS I A 

Xso/chord  =  .»  Yea/chord  =  ■  738 


14.  t 


INCHES 


WAKE  .  SUR.V E Y  DATA 

TEST  ill  RUN  UL 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  c  MU 


PRESSURE  ~  PSIA 

Xas/chord  —  Yes/chord  -  •  7^7 

BETA  =  SO 


WAKE  SURVEY  DATA 
TEST  zer  RUN  £o_ 

ANGLE  OF  ATTACK  O  NOMINAL  JET  ANGLE  3 fL  C  MU 


14.4  14. S  14.C 


15. •  15.1 


PRESSURE  ~  PS I A 


Xss/chord  —  .7&3  Yes/chord  =  .67 3 


•17.5  -15.5  -11.5 


Figure  III-9. 


ALPHA  ~  DEGREES 


INCHES  Z  -  INCHES 


WAKE  SURVEY  DATA 

TEST  lar  RUN  £JL 

ANGLE  OF  ATTACK  _2_  NOMINAL  JET  ANGLE  c  MU  £0**6 


BETA  =  JJL 


t 

N 


Figure  Ill-lO. 


ALPHA  ~  DEGREES 


191 


WAKE  SURVEY  DA.IA 

TEST  2M1  RUN 

ANGLE  OF  ATTACK  _JL_  NOMINAL  JET  ANGLE  J£_  C  MU  «//JV 


PRESSURE 

Xsa/chord  —  *  Z*i7  Yes/ chord  =  » 1 


INCHES 


WAKE  SURVEY  DATA 

TEST  ££Z  RUN  5~3 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  2L2_  C  MU 


PRESSURE  ~  PS I A 

Xss/chord  =  <  37  Yea/chord  =  •  V  9  2 


r* 


i 


16 


WAKE  SURVEY  DATA 

TEST  ISI  RUN  JJk_ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  C  MU  *< '02 


PRESSURE 

Xse/chord  =  ft.  /Q.&  Yes/chord  =  .3 


WAKE  SURVEY  DATA 

TEST  LSI  RUN  ■ 

ANGLE  OF  ATTACK  _JL  NOMINAL  JET  ANGLE  JL Q.  C  MU  8.660 6 


INCHES  |  Z  ~  INCHES 


WAKE  SURVEY  DATA 
TEST  ZSJ  RUN  &L.  ' 

ANGLE  OF  ATTACK  _SL  NOMINAL  JET  ANGLE  C  MU  gj] 


Xaa/chord  **  1?/$  Yaa/chord  “ 

BETA  =  'O 


Figure  III-20. 

201 


ALPHA  ~  DEGREES 


ANGLE  OF  ATTACK 


WAKE  SURVEY  DATA 
TEST  331  RUN  _62_  - 
a  NOMINAL  JET.  ANGLE  _2£_  C  MU  'h_Z7G>a 


WAKE  SURVEY  DATA 

TEST  i£Z  RUN  _£1 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  _2 0_.  C  MU  /6.8VS9 


PRESSURE  "  PS I A 

Xee/chord  **  .  Yaa/chord  *  6 

BETA  *  /? 


Figure  III-22. 


ALPHA  ~  DEGREES 


INCHES 


WAKE  SURVEY  DATA 

TEST  2&2  RUN  i/  • 

ANGLE  OF  ATTACK  _2_  NOMINAL  JET  ANGLE  _3£.  C  MU 


Xse/chord  =  <  Yae/chord  =  - 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 
TEST  HI  RUN  AL  ' 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET.  ANGLE  j£_  C  MU 

ta.c 

1M 

?.* 


1.1 

t.i 

1.1 

11.*  14.*  14.S  IS.*  IS.*  lt.l  11.1  IT.* 

PRESSURE  -  PS I A 

Xae/chord  =  •  Zi>fS  Yee/chord  -  >  63  ^ 

BETA  » 


INCHES  2  ~  INCHES 


WAKE  SURVEY  DA.IA 
TEST  ZS1  RUN  Je±  ' 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  AQ_  C  M'J 


Xee/chord  =  •/ ^  Ysa/chord  =  .SW 
BETA  =  / O 


Figure  III-28. 
209 


ALPHA  ~  DEGREES 


INCHES  2  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  Z&l  RUN  I£_  ■ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  JO_  C  MU  M 


BETA  =  *Q 


210 


ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  _tQ_  C  MU  .99*77 


PRESSURE  ~  PSIA 

Xss/chord  =  •  368  Yss/chord  =  » S' 

beta  =  zo 


INCHES  2  ~  INCHES 


Figure  III-32. 


ALPHA  ~  DEGREES 


WAKE. SURVEY  . DATA 

TEST  2£Z  RUN  38 


cn 

UJ 

X 

u 


X 

N 


CO 

LU 

X 

u 


X 

M 


Figure  III-33. 


ALPHA  ~  DEGREES 


1 


214 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  MZ  RUN  J±l_ 

ANGLE  OF  ATTACK  _j&_  NOMINAL  JET  ANGLE  5£_  C  MU  Ag// 

13.5 


!*.• 


H.3S  14.4*  14.45  14.5*  14.55  14.6*  14.65  14.7* 


PRESSURE  ~  PSIA 

Xse/chor'd  =  ,B37  Yes/chord  =  .3?y 
BETA  =  20 


n 

N 


Figure  III-34. 


ALPHA  ~  DEGREES 


215 


217 


INCHES  2  "  INCHES 


WAKE.SURVEY  DATA. 

TEST  IL&I  RUN  _1 0_ 

ANGLE  OF  ATTACK  „0_  NOMINAL  JET  ANGLE  JJl  C  MU  L °&2 


BETA  »  3Q 


Figure  III-38. 
219 


ALPHA  ~  DEGREES 


WAKE. SURVEY. DATA 

TEST  2B7  RUN  _3.3 

ANGLE  OF  ATTACK  _fl_  NOMINAL  JET  ANGLE  _7£L  C  MU 


WAKE  SURVEY  DATA 

TEST  tBl  RUN  _L 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET  ANGLE  C  MU 


14.2  14.2  14.4  14.2  14.2  It. 2  It. 2  It. 4 


PRESSURE  ~  PS I A 


X»»/chord  =  /6733  Y»«/chord  =  »  ~722-2- 
BETA  =  30 


WAKE  SURVEY. DATA 

TEST  2£7  RUN 

ANGLE  OF  ATTACK  _J2_  NOMINAL  JET  ANGLE  ±0_  C  MU  d367 


BETA  «  3o 


Figure  m-45.  ^PHA  ~  DEGREES 

226 


WAKE  SURVEY  DATA 

TEST  2JLZ.  RUN  3 


X««/chord  =  $S?S  Y»«/chord 


T 


-M 


T 


1« 


t — r 


~i — i — i — r 


»  F 


-i — i — i — r 


i* 


zt 


Figure  III-48. 


ALPHA  ~ 


229 


WAKE  SURVEY  DATA 

TEST  2R1  RUN 


ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  90 

«*  "1 - 1 - 1 - - - I - — - T - 


C  MU  /-  88/6 


INCHES  2  ~  INCHES 


WAKE  SURVEY  DATA 
TEST  132  RUN  7 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  C  MU  L922>^ 


BETA  -  JUL. 


Figure  III-50. 


ALPHA  ~  DEGREES 


M 


WAKE  SURVEY  DATA 

TEST  ££/  RUN  JL 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  C  MU  L M V? 


WAKE  SURVEY  DATA 

TEST  ig7  run 

*  *  r-  - 

ANGLE  GF  ATTACK  JL-  NOMINAL  JET  ANGLE  C  MU  /.  9VdS 


INCHES  2  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  Z33  RUN  / O 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  JO.  C  .MU  A 8989 


14.925  14. *75  14.125  14.175 

PRESSURE  ~  PS I A 

X««/chord  *  Ye®/chord  **  •*>  7 


9  '  19  29 


ALPHA  ~  DEGREES 


Figure  III-53. 


WAKE  SURVEY  DATA 

TEST  2&1  RUN  n 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  lQ_  C  MU  3J3/9 


INCHES  Z  ~  INCHES 


WAKE. SURVEY  DATA 
TEST  2jg.7_RUN  JJl _ ' 

ANGLE  OF  ATTACK  NOMINAL  JET  '  ANGLES.  C  MU  X7VI 


PRESSURE  ~  PS I A 


X®e/chord  =  Yss/ehord  -  A  ^ £  ^6 

BETA  = 


ALPHA  ~  DEGREES 

Figure  III-56. 
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WAKE  SURVEY  DATA 


INCHES  Z  ~  INCHES 


WAKE  SURVEY  DATA 

TEST  RUN  JJL 


■  Xaa/ chord  -  V  aa/ chord  - 

BETA  -  vr 


239 


INCHES 


WAKE  SURVEY  DATA 
TEST  2&Z  RUN  /6 

ANGLE  OF  ATTACK  _Q_  NOMINAL  JET  ANGLE  fJL  C  MU  7.yj/ 


13.75  14. ••  14.25  14.5*  14.75  IS. «•  1S.2S  15.5* 


PRESSURE  ~  PS I A 

X®»/  chord  55  ^  8 S' 7  Ye®/ohord  “  ‘382 
BETA  «  70 


WAKE  SURVEY  DATA 
TEST  231  R'JM  JJL 


V/A*<E  SURVEY  DATA 
TEST  ZB7  RUN  JA. 


INCHES 


WAKE  SURVEY  DATA 
TEST  232  RUN  JUL 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  C  MU  7. 7 8 S' 9 


BETA  «  To_ 


243 


INCHES  Z  "  INCHES 


WAKE  SURVEY  DATA 

TEST  £81  RUN  _2£ 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  JO_  C  MU  7. (,79 9 


X«»/ohord  “  2.1  S 6  Y©*/ehord  =  *  9*/ 

SETA  =  70 


244 


INCHES 


WAKE  SURVEY  DA  LA. 

TEST  231  RUN  3L 

ANGLE  GF  ATTACK  _P_.  NOMINAL  JET  ANGLE  9&_  C  MU  ".y*/3 


X««/ chord  —  /■  8 S~7  Y*«/ chord  ~  » 3  8  2. 


BETA  -  70 


Figure  III-64. 
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ALPHA  ~  DEGREES 


INCHES 


r  . .  ■ 

WAKE  SURVEY  DATA 

TEST  liLZ  RUN  LL. 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET  .ANGLE  1Q_  C  m /L*£89 


hi 

!HI 

— - 

i 

! 

mi 

h 

14.2  14.4  14.2  14.2  IS. 2  IS.  2  IS. 4 

PRESSURE  ~  PS I A 

Xaa/chord  **  X 899  Yaa/ehord  “  •  ZbiT 
BETA  =  70 


246 


WAKE  SURVEY  DATA, 

TEST  2MI  RUN  23 

ANGLE  OF  ATTACK  _0_  NOMINAL  JET  ANGLE  iSL.  C  MU  "  200/ 


Figure  III-67. 


ALPHA  ~  DEGREES 


248 


INCHES  2  ~  INCHES 


WAKE.  SURVEY  DATA 

TEST  23J  SUN  J2JT 

ANGLE  OF  ATTACK  _2_.  NOMINAL  JET  ANGLE  J&L  C  MU  //.vo93 


Xae/ chord  »  2,  idT 6  Yaa/  chord  *  J*.  V  V 

BETA  «*  70 


Figure  III-68.  ALPHA  ~  DEGREES 


249 


;KE  SURVEY  DATA 

ESI  2&Z  RUN  2JL 

CL  NOMINAL  JET  ANGLE  £Q.  Z  M'J  Lt&SZ 


Xea/chord 


PRESSURE  ~  PS I A 

**  ^  7  Y«»/chord  15  « 

BETA  =  70 


mm 


M 


i 


INCHES  2  "  INCHES 


WAKE  SURVEY  DATA 

TEST  287  RUN 

ANGLE  OF  ATTACK  NOMINAL  JET  ANGLE  _90  C  MU  <V.y7/7 


PRESSURE  "  PS I A 

Xaa/chord  “  A  Yea/  chord  “  •  2  6  £ 

BETA  =  ..?<>■ 


Figure  III-70. 


ALPHA  "  DEGREES 


■.KE  SURVEY  DATA 

ESI  ?31  RUN  _£!. 

o.  NOMINAL  JET  ANGLE  JLQ_  C  MU  8JI-.2.0  98 


DATE 

ILMED 


